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Abstract  of  Dissertation  Presented  to  the  Graduate  School 
of  the  University  of  Florida  in  Partial  Fulfillment  of  the 
Requirements  for  the  Degree  of  Doctor  of  Philosophy 

STRAIN  RELIEF  MODEL  FOR  THE  InxGai-xAs/GaAs  SYSTEM  AND  ITS 
APPLICATION  TOWARD  THE  GROWTH  OF  LOW  DISLOCATION 

DENSITY  InxGa  i -xAs 

By 

Viswanath  Krishnamoorthy 
August  1992 

Chairman:  Dr.  Robert.  M.  Park 

Major  Department:  Materials  Science  and  Engineering 
Strain  relaxation  occurring  in  the  InxGai-xAs/GaAs 
heteroepitaxial  system  as  a function  of  InxGai-xAs  composition  was 
studied  extensively  using  transmission  electron  microscopy,  to 
observe  dislocation  evolution,  and  high  resolution  x-ray  diffraction,  to 
determine  the  extent  of  strain  relaxation.  Both  single  InxGai-xAs/ 
GaAs  heterostructures  and  stacked  InxGai-xAs  multiple  layer  systems 
were  analysed. 

In  the  case  of  stacked  multiple  layer  InxGai-xAs/GaAs  systems, 
it  was  discovered  that  low-dislocation  density  InxGai-xAs  top  layers 

(layers  on  top  of  the  stacked  structures)  can  be  grown  provided  that 
the  compositional  difference  between  two  adjacent  InxGai-xAs  layers 
in  the  stack  is  maintained  below  a determined  critical  compositional 
difference  of  Ax=0.18,  taking  partial  lattice-relaxation  properly  into 
account.  For  this  case  it  was  determined  that  dislocation  evolution  is 
confined  to  the  GaAs  substrate  and  the  InxGai-xAs  layers  underlying 
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the  top  InxGa  i -xAs  layer  in  the  stack,  the  top  InxGa  i -xAs  layer  being 
essentially  dislocation-free.  This  phenomenon  is  attributed  to  a 
monotonic  increase  in  the  yield  strength  of  InxGai-xAs  at  the 
appropriate  growth  temperatures  with  increasing  values  of  x.  Such 
behaviour  appears  to  persist  up  to  an  InxGa  i -xAs  composition  of 
approximately  x=0.5,  whereupon  a further  increase  in  composition 
(increase  in  x value)  results  in  dislocation  evolution  in  the  top  layer  of 
the  stack.  It  is  postulated  that  the  yield  strength  of  InxGai-xAs 
decreases  at  elevated  temperatures  (growth  temperatures)  with 
increasing  values  of  x beyond  x=0.5. 

It  was  determined  that  extremely  low-dislocation  density 
InxGai-xAs  material  can  be  grown  on  GaAs  using  the  stacked  structure 
approach  as  evidenced  by  etch  pit  analysis.  For  example,  dislocation 
densities  of  1-2x10)3/ cm2  ancj  5-6x1  O^/cm^  were  recorded  from 
In0.35Gao.65As  and  Ino.48Gao.52As  top  layers,  respectively.  Such 
In^ai-xAs  alloys  would  be  potentially  suitable  for  the  fabrication  of 
photonic  devices  operating  at  1.3pm  (x=0.35)  and  1.55pm  (x=0.48). 
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CHAPTER  1 


INTRODUCTION 
1.1  Motivation  and  Objectives 

The  incorporation  of  semiconductor  heterostructures  into 
electronic  and  photonic  devices  has  in  recent  years  led  to  dramatic 
improvements  in  device  operating  characteristics.  For  instance, 
heteroemitter  bipolar  transistors  display  higher  current  gain  than 
conventional  homojunction  bipolar  transistors,  modulation-doped 
field -effect  transistors  or  high-electron  mobility  transistors  (HEMTs) 
operate  at  higher  speeds  than  do  conventional  field -effect  transistors 
and  double-heterostructure  diode  lasers  have  greatly  reduced 
threshold  current  densities  compared  to  their  homojunction 
counterparts. 

High  performance  electronic  and  photonic  devices  that 
incorporate  heterostructures,  such  as  those  mentioned  above,  have 
been  realized  due  to  recent  advances  in  epitaxial  growth  techniques 
such  as  molecular  beam  epitaxy  (MBE)  and  metal-organic  chemical 
vapor  deposition  (MOCVD).  These  advanced  epitaxial  growth 
techniques  are  low-temperature  in  nature  and  permit  atomic  level 
control  over  the  deposition  process  enabling  heterostructures  to  be 
grown  having  extremely  abrupt  interfaces. 
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A particular  advantage  of  modem  epitaxial  growth  techniques  is 
that  they  can  be  employed  to  integrate  materials  having  different 
properties  in  a monolithic  fashion  within  limits  imposed  by  factors 
such  as  lattice-mismatch  and  thermal  stability.  An  example  of 
materials  integration  would  be  the  GaAs/Si  epitaxial  system,  where 
photonic  devices  fabricated  in  GaAs  could  be  integrated  monolithically 
with  Si  electronic  devices.  A second  example  of  materials  integration 
concerns  the  subject  of  this  dissertation,  namely  the  InxGai-xAs/GaAs 
epitaxial  system. 

The  work  to  be  described  was  motivated  by  the  idea  of 
monolithically  integrating  photonic  devices  fabricated  in  InxGai-xAs 
with  GaAs  integrated  circuits  (ICs).  Such  a scheme  is  interesting  in 
principle  due  to  the  following  reasons:  Firstly,  GaAs  IC  technology  is 
rapidly  advancing  to  the  point,  for  example,  that  high-speed  GaAs  logic 
elements  have  been  incorporated  in  Cray  computers.  And  secondly, 
the  InjcGai-xAs  ternary  system  has  a direct-bandgap  range  from 
1.42eV  (x=0)  to  0.36eV  (x=l).  This  bandgap  range  would  correspond 
to  photonic  devices  operating  in  the,  0.87pm  (x=0)  to  3.4pm  (x=l), 
wavelength  range.  Such  a wavelength  regime  encompasses  several 
important  wavelengths.  For  instance,  wavelengths  of  1.3  and  1.55pm 
represent  the  dispersion  and  absorption  minima,  respectively,  in  the 
case  of  silica  glass  fibers  and  also  fluoride -based  fibers  appear  to  have 
an  attenuation  minimum  around  2pm. 

Hindering  the  growth  of  device  quality  InxGai-xAs  alloys  on 
GaAs,  however,  are  the  large  lattice-mismatch  strains  associated  with 
this  material  system  which  lead  to  dislocation  formation  and  thus  the 
degradation  of  photonic  devices.  Control  over  dislocation  propagation. 
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consequently,  is  of  paramount  importance  and,  hence,  fundamental 
studies  concerning  this  subject  are  required  such  that  the 
InxGai-xAs/  GaAs  material  integration  can  become  a reality. 

The  presently  described  work  had  two  principal  objectives. 
Firstly,  the  mechanisms  responsible  for  strain  relief  in  the 
InxGai-xAs/GaAs  heteroepitaxial  system  were  to  be  identified  as  a 
function  of  InxGai-xAs  composition,  covering  as  wide  a composition 
range  as  possible.  Secondly,  a composition  range  was  to  be  identified 
over  which  low-dislocation  density  InxGai-xAs  material  could  be 
grown  either  directly  or  indirectly  on  GaAs. 

Prior  to  outlining  the  scope  of  the  present  work,  a discussion 
section  follows  on  the  literature  as  it  pertains  to  this  dissertation. 

1.2  Background 

Lattice-mismatch  strain  presents  the  greatest  obstacle  towards 
the  growth  of  high  quality  InxGai-xAs  on  GaAs.  The  lattice  parameter 
of  InxGai-xAs  is  larger  than  that  of  GaAs  with  the  mismatch  strain 
reaching  7%  when  x=l  [1],  Therefore,  in  the  InxGai-xAs/GaAs 
system,  the  epilayer  experiences  in-plane  biaxial  compression  since 
the  epilayer  lattice  tries  to  conform  to  that  of  the  substrate  until  an 
elastic  limit  is  reached  and  strain  relaxation  occurs.  The  following 
section  reviews  the  published  literature  on  strain,  strain  relief  and 
dislocation  evolution  in  lattice -mismatched  systems  in  general,  with 
particular  reference  to  the  InxGai-xAs/GaAs  material  system.  Also, 
various  buffer  structures  used  to  block  dislocations  in  such  lattice- 
mismatched  systems  are  reviewed. 
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1.2.1  Critical  Thickness  in  Lattice-Mismatched  Systems 

In  the  case  of  heteroepitaxial  systems  in  which  the  lattice- 
mismatch  is  small,  the  first  atomic  layers  that  are  deposited  during 
epitaxial  growth  will  be  strained  to  match  the  substrate  lattice 
(pseudomorphic  growth)  and  a coherent  interface  is  formed. 

However,  as  the  thickness  of  the  epilayer  increases,  the  homogeneous 
strain  energy  in  the  epilayer  becomes  large  enough  to  allow  the 
nucleation  of  misfit  dislocations  at  the  heterointerface.  The  thickness 
of  the  epilayer  at  which  misfit  dislocations  are  nucleated  corresponds 
to  the  critical  thickness  for  that  epilayer/substrate  material  system.  In 
this  section  several  models  used  to  predict  critical  thickness  are 
discussed,  and  published  experimental  data  on  critical  thicknesses  in 
the  InxGai-xAs/GaAs  material  system  are  presented. 

Van  der  Merwe  [21  proposed  a model  to  predict  the  critical 
thickness  of  a lattice-mismatched  epilayer/  substrate  system  by 
equating  the  areal  strain  energy  density  in  the  epilayer  to  the 
interfacial  energy  per  unit  length  between  the  epilayer  and  the 
substrate  arising  as  a result  of  the  generation  of  misfit  dislocations. 

The  interfacial  energy  was  evaluated  using  the  Peierls-Nabarro 
potential  to  represent  the  interatomic  forces  at  the  interface.  Such  a 
calculation  of  the  interfacial  energy,  however,  becomes  significantly 
inaccurate  when  the  misfit  dislocation  spacing  is  greater  than  twice 
the  epilayer  thickness  leading  to  a conservative  estimate  of  the  critical 
thickness. 
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Using  a different  approach  than  that  employed  by  Van  der 
Merwe  [2],  Matthews  [3-5]  calculated  the  critical  thickness  by 
minimizing  the  elastic  (coherent)  and  plastic  (dislocation)  energy 
terms.  This  is  equivalent  to  a mechanical  equilibrium  of  forces  acting 
on  a dislocation.  Matthews  [3-5]  assumed  that  the  lattice -mismatch  is 
accommodated  by  a regular  array  of  non-interacting  edge  dislocations 
and  that  a misfit  dislocation  is  formed  by  the  bowing  of  grown-in 
threading  dislocations,  such  that  a segment  of  the  dislocation  line  lies 
parallel  to  the  interface.  The  forces  acting  on  a threading  dislocation 
are  the  line  tension  and  the  force  due  to  mismatch  stress  at  the 
interface.  In  the  initial  coherent  layer  the  line  tension  force  is  larger 
than  the  force  due  to  the  mismatch  stress.  As  the  epilayer  thickness 
increases,  the  force  due  to  the  mismatch  stress  increases  and  at  the 
critical  thickness  is  equal  to  the  line  tension  force.  This  mechanism, 
however,  becomes  inaccurate  for  small  lattice-mismatch  strains  due  to 
the  low  density  of  grown-in  threading  dislocations. 

Recently,  People  and  Bean  [6,7]  followed  a procedure  similar  to 
that  of  Van  der  Merwe  [2]  to  evaluate  critical  thicknesses  of  lattice- 
mismatched  systems.  These  authors  assume  that  the  growing  epilayer 
is  initially  free  of  threading  dislocations  and  that  the  interfacial  misfit 
dislocations  are  generated  when  the  areal  strain  energy  density  in  the 
epilayer  exceeds  the  self  energy  of  an  isolated  dislocation  at  the 
interface.  The  energy  balance  is  considered  over  a region  of  the 
coherently  strained  lattice  equal  to  the  width  of  the  dislocation  strain 
field.  There  are  three  problems  associated  with  People  and  Bean’s 
[6,7]  calculations.  Firstly,  the  energy  of  an  isolated  screw  dislocation 
in  an  infinite  crystal  lattice  is  used  to  represent  the  energy  of  a 


6 


dislocation  in  a finite  epilayer.  Secondly,  the  interface  consisting  of 
two  vastly  different  materials  (in  terms  of  mechanical  properties)  is 
considered  to  have  similar  properties  as  that  of  a bulk  epilayer 
material.  Thirdly,  an  effective  interfacial  width  of  only  four  to  five 
<110>  atomic  spacings  is  used  to  determine  the  energy  of  the 
dislocation  which  is  inaccurate  for  small  lattice-mismatched  systems. 
The  aforementioned  factors  lead  to  a gross  overestimation  of  the 
critical  thickness. 

Several  experimental  studies  concerning  critical  thickness 
measurements  in  the  InxGai-xAs/GaAs  system  performed  using 
various  analytical  techniques  have  been  published  in  the  literature  [8- 
13]  and  the  reported  data  are  illustrated  in  Figure  1-1  together  with 
the  theoretical  predictions  of  Van  der  Merwe  [2],  Matthews  [3-5]  and 
People  and  Bean  [6,7].  Orders  and  Usher  [9]  determined  lattice 
parameters  of  InxGai-xAs  (x=0.07,0. 14  and  0.25)  layers  of  various 
thicknesses  grown  on  GaAs  substrates  using  x-ray  rocking  curve 
analysis  and  denoted  the  critical  thickness  (hc)  to  be  the  layer 
thickness  at  which  the  lattice  parameter  changed  abruptly  to  its 
relaxed  value.  Their  [9]  determination  of  critical  thicknesses  seemed 
to  agree  with  the  theoretical  predictions  of  People  and  Bean  [6,7] 
(Figure  1-1)  for  InxGai-xAs  compositions  of  x=0.14  and  x=0.25. 

Another  critical  thickness  study  was  reported  by  Fritz  et  al.  [10] 
who  performed  photoluminescence  (PL)  microscopy  on  Ino.2Gao.8As 
layers  of  various  thicknesses  grown  on  GaAs.  The  PL  micrographs 
showed  the  presence  of  dark-line  defects  for  epilayer  thicknesses 
greater  than  20nm.  As  can  be  seen  from  Figure  1-1  this  value  of 
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Figure  1-1.  InxGai-xAs/GaAs  epilayer  critical  thickness  data  reported 
by  various  authors  [8-13]  as  a function  of  InxGai-xAs 
composition.  Also  shown  are  theoretical  plots  based  on 
three  models  [2-7]. 
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critical  thickness  lies  between  the  curves  predicted  by  Matthews  [3-5] 
and  People  and  Bean's  [6,7]. 

Nakao  and  Yao  [8]  performed  in-situ  reflection  high  energy 
electron  diffraction  (RHEED)  during  molecular  beam  epitaxial  (MBE) 
growth  of  InxGai-xAs  layers  (x=0.23,  0.66,  0.78  and  1)  on  GaAs  and 
measured  the  surface  lattice  parameter  of  the  epilayers  from  the 
spacing  between  reciprocal  rods  as  a function  of  epilayer  thickness.  At 
the  critical  thickness,  the  surface  lattice  parameter  abruptly  changes 
from  the  substrate  (coherent)  to  the  bulk  epilayer  (relaxed)  value. 
Similar  investigations  were  conducted  by  Price  [11]  for  InxGai-xAs 
layer  compositions  of  x=0.3,  0.38  and  0.5  and  by  Berger  et  al.  [12]  for 
x=0.35.  All  the  experimental  values  are  shown  in  Figure  1-1  and  the 
data  points  lie  between  Matthew’s  [3-5]  and  Van  der  Merwe's  [2] 
theoretical  models  for  the  evaluation  of  critical  thickness.  In  a recent 
study.  Price  [14]  showed  that  the  critical  thickness  of  InxGai-xAs  is  a 
strong  function  of  the  temperature  employed  in  epitaxial  growth. 

Price  [14]  observed  that  for  an  Ino.35Gao.65As  epilayer  grown  on  GaAs 
by  the  MBE  technique,  the  critical  thickness  increases  with 
decreasing  growth  temperature.  Price  [14]  also  showed  that  the 
surface  lattice  parameter  shows  a gradual  increase  with  increasing 
thickness  until  a critical  growth  temperature  is  reached  when  the 
lattice  parameter  increases  abruptly  at  a critical  thickness.  Ideally, 
during  the  epitaxial  growth  process,  the  epilayer  grows  one  atomic 
layer  at  a time  (2D  growth).  In  order  to  achieve  this  ideal  growth 
condition,  the  surface  mobility  of  the  atoms  should  be  large  so  that 
they  can  migrate  to  an  atomic  ledge  on  the  surface  and  allow  for 
uniform  coverage  throughout  the  monolayer.  However,  if  the  lattice- 
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mismatch  is  large  (>2%)  [11],  the  mismatch  strain  might  significantly 
decrease  the  surface  mobility.  Also,  if  the  growth  temperature  is  low, 
the  surface  mobility  of  atoms  is  significantly  reduced.  Therefore  in 
scenarios  of  large  lattice  mismatch  strains  and/or  low  growth 
temperatures,  island  (3D)  growth  occurs.  Island  growth  allows  for  a 
larger  areal  strain  energy  than  2D  growth  and  hence  the  observed 
larger  critical  thickness  at  lower  growth  temperatures.  The  growth 
temperature  at  which  the  surface  lattice  parameter  shows  an  abrupt 
change  is  denoted  as  the  3D-2D  growth  transition  temperature.  Other 
published  results  [12,15-17]  of  similar  RHEED  studies  on  InxGai-xAs 
epilayers  as  a function  of  the  growth  temperature  report  similar 
observations  and  use  the  2D-3D  growth  transition  temperature 
concept  to  explain  the  observations. 

In  summary,  the  experimental  data  reported  in  the  literature 
suggest  that  the  onset  of  strain  relief  (critical  thickness)  in  the 
InxGai-xAs/GaAs  system  does  not  conform  precisely  to  any  of  the 
available  theoretical  models  concerning  the  prediction  of  critical 
thickness.  Also,  the  critical  thickness  of  InxGai-xAs/GaAs  epilayers 
appears  to  be  sensitive  to  the  growth  temperatures  employed. 

1.2.2  Strain  Relief  Mechanisms 


Strain  relief  occurs  in  lattice -mismatched  systems  by  the 
evolution  of  dislocations  when  the  epilayer  thickness  reaches  a critical 
value.  In  this  section,  the  published  literature  on  proposed  models  for 
strain  relaxation  and  experimental  results  concerning  the  InxGai-xAs/ 
GaAs  material  system  are  discussed. 
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Matthews  and  Blakeslee  [18]  proposed  a model  to  predict  the 
extent  of  strain  relaxation  by  considering  an  equilibrium  of  stresses 
acting  on  a threading  dislocation.  The  stress  exerted  by  the  lattice 
mismatch  strain  on  a dislocation  (oe)  in  its  glide  plane  and  the 
effective  stress  due  to  the  dislocation  line  tension  (ai)  are  given  by 
2 pe  (1  + v) 


and 

1 2n  1 - v h/b 

where  p is  the  shear  modulus,  v is  the  Poissons  ratio,  b is  the  Burgers 
vector  and  h is  the  thickness  of  the  epilayer.  The  two  stresses  are 
equated  to  give  the  residual  strain  in  the  epilayer 

e b I1'?)  1 + ln(f) 

Ceq  4k  (1  + v)  h 

where  eeq  is  the  equilibrium  residual  strain  in  the  epilayer. 

A second  model  proposed  by  Tsao  et  al.  [19]  suggests  that  the 
driving  stress  for  dislocation  motion  is  the  difference  between  the 
stress  due  to  misfit  strain  (oe)  and  the  effective  stress  due  to  the 
dislocation  line  tension  (aj)  and  an  observable  strain  release  should 
occur  only  if  this  "excess"  stress  (aex)  exceeds  a critical  value  which 
depends  on  the  growth  temperature. 


<*ex  = Oe  - ai 


1 1 


Using  the  excess  stress  a shift  from  the  equilibrium  residual  strain 
value  was  given  by 

a^d  -v) 

em  " ^ + 2|l(l  + V) 

Since  the  value  of  Cex  is  unknown,  a number  of  plots  of  em  versus  h 
for  various  values  of  the  excess  stress  are  used  to  obtain  a fit  to 
experimental  data  unique  to  a lattice-mismatched  system. 

Another  model  for  strain  release  was  proposed  by  Dodson  and 
Tsao  [20,21]  assuming  that  strain  relaxation  is  limited  by  the  kinetics 
of  dislocation  movement  and  multiplication.  Without  going  into  the 
mathematics,  the  final  result  is  a differential  equation  given  by 

^ = cn2  [f-  Y(t)  - eeq]2  «t)  + y0] 

where  y(t)  is  the  time-dependent  strain  release,  eeq  is  the  equilibrium 
strain,  yo  is  a source  term  which  accounts  for  the  density  of  threading 
dislocations  in  the  substrate,  f is  the  misfit  strain  and  C is  a 
phenomenological  parameter  which  contains  the  dislocation-glide 
velocity.  Due  to  the  term,  C,  the  strain  relaxation  should  depend  on 
the  growth  velocity,  i.e.  the  time  the  specimen  is  held  at  the  growth 
temperature. 

Using  a different  approach,  Maree  et  al.  [22]  explained  strain 
relaxation  in  the  epilayer  by  means  of  surface  nucleation  and  expansion 
of  dissociated  half  loops  while  taking  into  account  the  misfit  stress 
field  acting  on  expanding  loops.  The  half  loop  energy  (Eh)  increases 
with  increasing  loop  radius  (r),  reaches  a maximum  value  when 
dEh/dr=0  at  a critical  radius  rc  and  then  decreases  with  further 
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dislocation  loop  expansion.  In  the  case  of  an  in-plane  compressive 
strain  (as  in  the  InxGai-xAs/GaAs  system),  the  residual  strain  in  the 
epilayer  (e)  is  given  by 


In  order  to  determine  the  model  that  best  predicts  the  extent  of 
strain  relaxation,  Drigo  et  al.  [23]  performed  a study  of  the  residual 
strain  in  MBE  grown  InxGai-xAs  (x=0.035,  0.08,  0.10  and  0.15) 
epilayers  of  various  thicknesses  grown  on  GaAs.  The  parallel  or  in- 
plane lattice -mismatch  residual  strain  was  determined  using 
Rutherford  backscattering/ channeling  spectrometry.  The  results 
showed  that  the  measured  residual  in-plane  strain,  represented  by 
data  points  and  a best  fit  to  the  data  denoted  by  a dotted  line  in  Figure 
l-2a,  is  much  larger  than  predicted  by  the  Matthews  and  Blakeslee 
[18]  equilibrium  model  represented  by  the  solid  line  in  Figure  l-2a. 

The  experimentally  determined  values  for  residual  strains  did  not 
agree  with  those  predicted  by  Tsao  et  al.  [19]  (Figure  l-2b).  The 
residual  strain  data,  however,  appears  to  fit  the  energy-balance  model 
proposed  by  Maree  et  al.  [22]  based  on  the  nucleation  of  half-loops 
(Figures  l-2b).  An  annealing  study  at  the  growth  temperature  was 
conducted  in  order  to  observe  if  further  strain  relief  occurred  with 
increasing  the  time  at  the  growth  temperature  as  proposed  by  Dodson 
and  Tsao  [20,21].  The  results  showed  no  further  strain  relief  thus 
discounting  Dodson  and  Tsao’s  model  [20,21]  for  strain  relief  in  the 
InjcGa  i -xAs/ GaAs  system.  Strain  relief  data  obtained  from  InxGai-xAs 
layers  by  Kamigaki  et  al.  [24]  for  x=0. 082-0.2  and  Orders  and  Usher  [9] 
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(a) 


Figure  1-2.  A plot  of  the  experimental  data  from  studies  [9,22-24]  on 
the  in-plane  residual  strain  (e)  as  a function  of  InxGai-xAs 
layer  thickness  (h)  and  its  correlation  to  proposed  models 
for  strain  relaxation. 

a)  The  dashed  line  represents  a best  fit  to  the 
experimental  data  while  the  solid  line  represents  the 
equilibrium  model  [18]  showing  that  the  model 
overestimates  the  extent  of  strain  relaxation;  b)  The  solid 
lines  represent  the  predictions  of  Tsao  et  al.'s  [19]  model 
for  various  values  of  the  excess  stress  and  show  no 
correlation  to  the  experimental  data.  The  dot-dashed  line 
represents  the  results  of  the  energy-balance  model  for  the 
nucleation  of  half-loops  [22]  and  provides  a good 
correlation  to  the  experimental  data  points.  (From  Drigo 
et  al.  [23]) 
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for  x=0.25  by  x-ray  diffraction  and  Maree  et  al.  [22]  for  x=0.07  by  ion 
blocking  are  also  shown  in  Figures  l-2a  and  l-2b.  These  data  points 
also  fit  the  energy  balance  model  based  on  the  nucleation  of  half  loops. 

In  summary,  as  far  as  the  literature  is  concerned  the  strain 
relaxation  model  based  on  half  loop  nucleation  and  expansion  appears 
to  best  represent  the  extent  of  strain  relaxation  in  the 
InxGai-xAs/GaAs  material  system. 

1.2.3  Dislocation  Evolution  in  the  In^Gai  -vAs/GaAs  Material  System 

In  order  to  nucleate  dislocations,  either  the  stress  in  the 
epilayer/ substrate  system  requires  to  reach  a value  that  is  favourable  to 
nucleate  dislocations  for  homogeneous  nucleation  or  there  has  be  a 
dislocation  source  such  as  a stress  concentrator  to  allow  for  the 
heterogeneous  nucleation  of  dislocations.  There  are  various  possible 
sources  that  could  act  as  stress  concentrators  for  dislocation  formation 
in  lattice-mismatched  systems. 

The  first  candidates  for  dislocation  sources  as  suggested  by 
Matthews  et  al.  [3]  are  threading  dislocations  that  replicate  from  the 
substrate  and  penetrate  the  strained  epilayer  (Figure  l-3a).  Such 
dislocations  will  produce  misfit  dislocations  when  Oex>0.  he.  the 
mismatch  stress  is  greater  than  the  dislocation  line  tension.  The 
relaxation  rate  is  proportional  to  the  velocity  of  the  threading 
dislocation  gliding  due  to  the  excess  stress. 
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Figure  1-3.  Schematic  diagram  of  possible  sources  of  misfit 
dislocations. 

a)  Generation  of  misfit  dislocations  from  grown-in 
threading  dislocations;  b)  Homogeneous  nucleation  of 
dislocation  half  loops  at  surface  steps;  c)  Heterogeneous 
nucleation  of  dislocations  at  stress  concentrators  such  as 
an  oxide  inclusion  which  could  act  as  dislocation  mills. 
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Another  possibility  is  the  homogeneous  nucleation  of  dislocation 
half  loops  at  the  surface  (Figure  l-3b).  In  order  to  activate  such  a 
source  at  the  surface,  the  physical  stress  at  the  surface  must  be  large 
enough  to  nucleate  dislocation  half  loops.  This  process  has  been 
addressed  by  several  authors  [25-27]  and  the  concensus  is  that  a 
lattice-mismatch  strain  of  several  percent  would  be  required  to 
nucleate  dislocations  at  such  surface  steps. 

A third  possibility  would  be  the  local  formation  of  dislocation 
mills,  similar,  in  principle,  to  a Frank-Read  source  [26]  (Figure  l-3c). 
The  mechanism  of  dislocation  mills  involves  the  expansion  of  a 
dislocation  structure  until  an  unbound  dislocation  loop  pinches  off, 
leaving  behind  the  original  dislocation  structure  which  is  then  free  to 
generate  more  dislocation  loops. 

Fitzgerald  et  al.  [28,29]  performed  an  extensive  study  of 
dislocation  generation  specifically  on  1pm  thick  InxGai-xAs  (x=0.12) 
epilayers  grown  on  GaAs  using  plan-view  and  cross-sectional  TEM 
analyses.  Plan-view  TEM  analysis  of  the  InjcGa  i -xAs /GaAs 
heterointerface  showed  a dislocation  network  made  up  of  60° 
a/2<110>  mixed  dislocations.  Since  InxGai-xAs  has  a zincblende-type 
structure,  perfect  dislocations  in  the  lattice  have  an  a/2<110>  type 
Burgers  vector.  Of  the  three  main  types,  namely,  edge,  screw  or  60° 
mixed  dislocations,  the  edge  type  has  the  highest  core  energy  and  the 
screw  type  cannot  relax  a tetragonal  lattice  mismatch,  leaving  the  60° 
mixed  type  as  the  most  probable  dislocation  to  form  in  InxGai-xAs 
layers  [30]. 

Furthermore,  cross-sectional  TEM  analysis  of  the 
heterostructure  showed  the  presence  of  dislocation  half-loops  in  the 
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GaAs  substrate  and  an  absence  of  dislocations  in  the  InxGai-xAs 
epilayer.  The  three  types  of  dislocation  sources  mentioned  earlier, 
namely,  threading  dislocations,  homogeneous  nucleation  at  the  surface 
steps  and  heterogeneous  nucleation  or  dislocation  mills  at  the  surface 
which  when  active  would  cause  threading  dislocation  formation  in  the 
epilayer  are  obviously  not  active  in  the  InxGai-xAs  (x=0. 12)/GaAs 
system  since  dislocations  are  not  observed  in  the  InxGai-xAs  epilayer. 
Fitzgerald  et  al.  [28,291  attributed  the  absence  of  dislocations  in  the 
epilayer  to  a larger  misfit  strain  in  the  epilayer  compared  to  the 
substrate.  However,  observations  in  other  lattice -mismatched  systems 
such  as  ZnSe/GaAs  [31]  and  GaAs/Si  [32],  where  the  misfit  strain  in 
the  epilayer  is  larger  than  that  in  the  substrate,  do  not  conform  to 
Fitzgerald  et  al.'s  [28]  explanation  since,  in  these  material  systems, 
dislocations  are,  in  fact,  observed  in  the  epilayer. 

Fitzgerald  et  al.  [28]  also  proposed  a model  to  predict  the  extent 
of  dislocation  propagation  into  the  substrate  by  considering  a balance 
of  forces  acting  on  a dislocation  in  the  GaAs  at  a distance  x along  {111} 
planes  from  the  heterointerface.  This  distance,  x,  is  correlated  to  a 
depth  of  penetration  of  a 60°  dislocation  into  the  substrate  by  straight 
forward  geometrical  considerations.  In  this  model,  a Peierls  force 
(Fp),  an  image  force  (Fi),  a line  tension  force  (Fi)  and  a force  from  the 
elastic  strain  in  the  substrate  (Fe)  act  on  a dislocation  to  prevent  it 
from  moving  further  into  the  substrate  and  repulsion  forces  (Ftf)  from 
60°  dislocations  having  the  same  Burgers  vector  act  to  move  the 
dislocation  into  the  substrate.  At  the  equilibrium  distance,  x, 

Fd  = Fe  + Fi  + Fj  + Fp. 
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Fitzgerald  et  al.  [28]  showed  that  the  calculated  value  for  penetration 
depth  of  a dislocation  into  the  substrate  agrees  with  the  experimental 
values  determined  using  XTEM. 

Chang  et  al.  [33,34]  conducted  PTEM  and  XTEM  analyses  on 
InxGai-xAs  (x=0.15,  0.3  and  0.5)  layers  grown  on  GaAs.  PTEM 
analysis  of  the  heterostructures  showed  that  the  dislocations  were 
mostly  60°  a/2<110>  mixed  type  in  character  at  the  Ino.15Gao.85As/ 
GaAs  heterointerface.  However,  analysis  of  the  Ino.3Gao.7As  and 
Ino.5Gao.5As  alloys  revealed  predominantly  edge  type  dislocations. 

The  increase  in  edge  type  dislocations  with  increasing  lattice- 
mismatch  was  attributed  to  island  growth  since  island  coalescence 
normally  leads  to  the  formation  of  edge  dislocations.  XTEM  analysis  of 
these  particular  InxGai-xAs/GaAs  structures  showed  the  absence  of 
dislocations  in  the  Ino.15Gao.85As  epilayer  and  the  presence  of  a 
large  number  of  threading  dislocations  in  the  Ino.3Gao.7As  and 
Ino.5Gao.5As  epilayers.  Chang  et  al.  [33]  suggested  a 2D-3D  growth 
mode  transition  due  to  large  lattice-mismatch  strains  (>2%)  to  explain 
the  appearance  of  threading  dislocations  in  x>0.15  epilayers.  Such  an 
explanation,  however,  does  not  appear  to  be  valid  since  studies  on 
material  systems  which  contain  much  smaller  lattice-mismatches  and 
therefore  have  a 2D  growth  mode,  for  example  the  ZnSe/GaAs  material 
system  [31],  show  threading  dislocations  in  the  epilayer  unlike 
InxGa  1 -xAs / GaAs  layers  with  similar  lattice-mismatch  strains.  Also, 
RHEED  and  TEM  studies  conducted  by  Tabuchi  et  al.  [35]  on  the 
InxGai-xAs/GaAs  system  (x=0.2,  0.3,  0.4,  0.5,  0.7  and  1)  show  that 
despite  an  observed  2D  growth  mode  for  InxGai-xAs  epilayers  with 
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x<0.5,  threading  dislocations  form  in  InxGai-xAs  epilayers  of  x>0.2. 
Therefore,  the  2D-3D  growth  mode  transition  mechanism  as  a means 
to  explain  the  onset  of  threading  dislocation  formation  in  InxGai-xAs/ 
GaAs  layers  appears  to  be  invalid  for  epilayers  with  x values  less  than 
0.5.  Tabuchi  et  al.  135]  stated  that  island  growth  occurs  in  InxGai- 
xAs/GaAs  systems  with  x>0.5.  Also,  Guha  et  al.  [36]  conducted  high 
resolution  XTEM  analysis  on  extremely  thin  layers  of  Ino.5Gao.5As 
grown  on  GaAs  and  directly  observed  the  formation  of  islands. 
Therefore,  it  appears  that  for  x>0.5,  island  formation  and  coalescence 
could  contribute  to  threading  dislocation  formation  in  InxGai-xAs 
epilayers. 

In  summary,  it  seems  clear  that  current  postulates  which  are 
designed  to  explain  the  nature  of  dislocation  sources  and  dislocation 
evolution  in  lattice-mismatched  systems  are  not  adequate  to 
completely  account  for  dislocation  evolution  in  the  InxGai-xAs/GaAs 
material  system. 

1.2.4  Common  Buffer  Structures  Used  to  Block  Dislocations 

Various  schemes  have  been  studied  by  various  authors  [37-40] 
with  a view  towards  reducing  the  threading  dislocation  density  in 
lattice-mismatched  InxGai-xAs  layers  grown  on  either  GaAs  or  InP. 

The  techniques  that  have  been  employed  to  minimize  the  effects  of 
lattice-mismatch  can  be  summarized  as  follows: 

a)  Growth  of  a continuously  graded  InvGai.YAs  buffer  layer  [38]  on  InP 
such  that  the  lattice  parameter  (corresponding  to  a layer  composition 
°f  In0.53Gao.47As)  adjacent  to  the  substrate  was  matched  to  that  of 
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InP  and  the  lattice  parameter  (and  composition)  at  the  top  of  the 
buffer  was  matched  to  that  of  an  Ino.77Gao.23As  epilayer  of  interest. 
Reduction  in  dislocation  density  via  use  of  the  buffer,  however,  was 
very  insignificant  (less  than  an  order  of  magnitude). 

b)  Step  graded  growth  where  thin  InxGai-xAs  layers  (-2000A  each) 
with  increasing  composition  (x=0.25,  x=0.5  and  x=0.75)  are  grown  on 
GaAs  in  a stacked  manner,  resulting  in  partial  confinement  of  the 
misfit  dislocations  at  intermediate  interfaces  [37,39].  Reductions  in 
dislocation  densities  by  up  to  2 to  3 orders  of  magnitude  were 
achieved  by  this  method  in  the  case  of  heteroepitaxy  of  InAs  on  GaAs. 

c)  Growth  on  patterned  GaAs  substrates  where  substrates  are  etched 
to  form  mesas  and  trenches  and  misfit  dislocations  form,  upon  growth 
of  an  Ino.05Gao.95As  layer,  by  the  nucleation  of  half-loops  which 
appear  to  glide  along  the  interface  leaving  a strain  relieving  dislocation 
behind  [41,42],  These  dislocations  glide  from  the  mesas  into  the 
trenches  but  appear  to  be  blocked  from  gliding  into  adjacent  mesas 
leaving  the  mesas  relatively  dislocation-free.  However,  InxGai-xAs 
layers  with  larger  lattice-mismatches  to  GaAs  require  to  be  studied  in 
order  to  justify  the  use  of  this  scheme  since,  as  has  been  detailed 
earlier,  such  small  x (x<0. 15)  InxGai-xAs  layers  can  be  grown  directly 
on  GaAs  without  the  introduction  of  threading  dislocations  in  the 
epilayer. 

d)  Growth  of  strained  layer  superlattices  (SLSs)  are  also  used  to  act  as 
barriers  to  dislocation  motion.  SLSs  are  the  most  widely  employed 
structures  designed  to  block  dislocations  and  hence  the  remainder  of 
this  section  is  devoted  to  the  principles  of  dislocation  blocking  by 
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SLSs,  in  general,  and  to  a discussion  on  the  experimental  studies 
specific  to  the  InxGai-xAs/GaAs  material  system. 

El-Masry  et  al.  [43]  have  studied  dislocation  interactions  within 
an  InxGai-xAs/GaAs  SLS  system  and  report  that  the  60°  mixed 
dislocations  are  bent  inside  the  SLS  and  their  direction  is  changed 
towards  one  of  the  <110>  directions  along  the  SLS  interfaces.  These 
authors  suggest  that  effective  bending  of  dislocations  requires  that  all 
threading  dislocations  experience  a gliding  force  known  as  the  Peach- 
Koehler  force  whose  magnitude  depends  on  the  strain  and  layer 
thickness. 

Matthews  and  Blakeslee  [18],  in  general,  have  suggested  that  a 
minimum  layer  thickness  in  the  SLS  is  required  in  order  to  bend 
threading  dislocations,  their  theory  being  based  on  the  existence  of  a 
finite  stress  in  individual  layers  of  the  SLS  acting  to  bend  over  and 
confine  threading  dislocations.  The  critical  minimum  thickness  is 
therefore  dependent  on  the  strain  between  the  individual  layers  in  the 
SLS. 

In  general  terms,  the  lattice-mismatch  in  a SLS  results  in  a 
biaxial  state  of  stress  for  the  individual  layers.  Assuming  isotropic 
behaviour,  the  shear  stresses  and  shear  strains  are  equal  to  zero  along 
the  {100}  planes.  However,  finite  resolved  shear  stresses  exist  in 
{111}  planes  along  <110>  directions  which  are  the  slip  planes  and  slip 
directions,  respectively,  for  60°  mixed  type  dislocations  which  are  the 
most  common  of  the  dislocation  types  observed  in  lattice-mismatched 
systems.  For  relatively  small  shear  stresses  dislocations  would 
continue  to  thread  into  the  epilayer,  hence  rendering  the  SLS 
ineffective;  however,  for  relatively  large  stresses  which  are  sufficient 
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to  bend  dislocations,  a significant  reduction  in  dislocation  density  can 
be  achieved  when  the  majority  of  the  dislocations  are  of  the  mixed 
type. 

Interactions  between  adjacent  dislocations  may  occur  if  the 
density  of  dislocations  penetrating  the  SLS  is  very  large  (lO^/cm^). 
These  interactions  depend  on  the  stress  field  around  the  dislocations 
and  the  distance  between  adjacent  dislocations.  If  an  attractive  stress 
field  exists  between  two  dislocations  they  can  interact  to  form  one 
dislocation  or  a half  loop,  thus  reducing  the  dislocation  density. 

Favourable  conditions  for  dislocation  reduction  would  be  realized 
when  1)  the  dislocation  is  bent  at  the  SLS  interface  and  propagates  to 
the  sample  edge,  2)  two  dislocations  interact  to  cancel  each  other  and 
form  a half  loop  and  3)  two  dislocations  react  to  form  a third 
dislocation  at  a node. 

Several  studies  [30,44,45]  of  dislocation  evolution  in 
InxGai-xAs/GaAs  (x=0.08-0.35)  SLSs  grown  on  GaAs  using  XTEM 
analysis  have  been  conducted  and  the  results  suggest  that  most  of  the 
dislocations  are  confined  to  the  SLS/GaAs  heterointerface  and  the 
GaAs  substrate,  with  only  a few  dislocations  penetrating  the 
superlattice.  Also,  it  appears  that  the  interactions  of  60°  misfit 
dislocations  lead  to  the  formation  of  edge  dislocations  at  the  nodes. 

The  use  of  such  low  In  content  (x<0.4)  InxGai-xAs  layers  in 
InxGa  i -xAs/ GaAs  superlattices  could  allow  for  the  growth  of  high 
quality  lattice -matched  layers  of  InxGai-xAs  with  x values  <0.2. 
However,  as  discussed  in  earlier  sections,  such  low  dislocation  density 
InxGai-xAs  layers  (x<0.2)  can  be  grown  directly  on  GaAs  without  the 
use  of  SLSs.  In  other  words,  SLS  buffer  layers  would  appear  to  be 
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unnecessary  for  dislocation  density  reduction  in  InxGai-xAs  layers  for 
x values  <0.2. 

1.2.5  Summary 

A detailed  review  of  published  literature  on  the  InxGai-xAs/GaAs 
system  suggests  that  currently  available  theoretical  models  for 
predicting  critical  thicknesses  in  lattice-mismatched  systems  are  not 
suitable  to  explain  the  onset  of  strain  relaxation  in  the  InxGai-xAs/ 
GaAs  material  system.  However,  the  extent  of  strain  relaxation  in  this 
system  appears  to  obey  the  theoretical  model  of  Maree  et  al.  [22]  based 
on  the  surface  nucleation  and  expansion  of  half  loops. 

Also,  published  studies  on  dislocation  evolution  in  the 
InxGai-xAs/GaAs  system  suggest  an  absence  of  dislocations  in  the 
epilayers  for  x<0.2.  However,  current  postulates  concerning  the 
nature  of  dislocation  sources  and  dislocation  evolution  are  not 
adequate  to  explain  dislocation  nucleation  and  evolution  in  this 
material  system. 


1.3  Scope  of  the  Present  Work 

The  present  work  is  detailed  in  five  following  chapters.  In 
chapter  2,  the  growth  techniques  employed  in  the  fabrication  of  single 
and  multilayer  InxGai-xAs  structures  on  GaAs  are  discussed  along  with 
the  analytical  techniques  used  to  characterize  the  various 
heterostructures. 
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Chapter  3 details  a dislocation  evolution  study  performed  on 
single  InxGai-xAs  layers  of  various  compositions  grown  on  GaAs. 
Dislocation  formation  is  explained  based  on  the  high  temperature 
mechanical  properties  of  the  epilayers  and  the  substrate.  A unique 
concept  of  "critical  composition"  is  proposed  based  on  a balance  of 
forces  model  which  includes  a surface  image  force  to  explain 
dislocation  evolution  in  the  InxGai-xAs/GaAs  material  system. 

In  chapter  4,  the  step  by  step  development  of  a multistage  strain 
relief  buffer  employing  conventional  dislocation  blockers,  i.e.  strained 
layer  superlattices,  for  the  growth  of  InxGai_xAs  layers  with  x values 
greater  than  0.2  is  discussed.  The  lattice-matching  between  individual 
layers  in  the  multistage  system  is  also  detailed  and  the  advantages  and 
disadvantages  of  using  such  a buffer  system  are  discussed. 

Chapter  5 details  a study  of  strain  relaxation  in  single 
InxGa  i -xAs/ GaAs  layers,  the  results  of  this  study  being  applied  towards 
the  development  of  a multilayer  system  which  permits  the  growth  of 
extremely  low-dislocation  density  InxGai-xAs  for  x values  up  to  0.5. 

The  multilayer  design  is  based  on  an  experimentally  determined 
critical  compositional  difference"  associated  with  adjacent  layers  of  a 
stacked  InxGai-xAs  multi-layer  system.  The  stacked  systems  consist 
of  relatively  thick  InxGaj-xAs  layers  and  do  not  employ  conventional 
dislocation  blockers,  i.e.,  strained -layer  superlattices. 

In  chapter  6,  a summary  of  the  important  conclusions  drawn 
from  this  study  is  provided  along  with  recommendations  for  future 
research  in  this  area. 


CHAPTER  2 

GROWTH  AND  CHARACTERIZATION  TECHNIQUES 

The  various  single  and  multiple  layer  InxGai-xAs/GaAs 
heterostructures  studied  in  this  work  were  grown  using  the  molecular 
beam  epitaxy  (MBE)  technique  by  Dr.  Y.W.  Lin  and  Bruce  Liu  of  our 
research  group,  while  the  heterostructures  were  designed  by  the 
author  and  characterized  using  various  analytical  techniques  including 
cross-sectional  and  plan-view  transmission  electron  microscopy 
(XTEM  and  PTEM),  high-resolution  x-ray  diffraction  (HRXRD), 
electron  microprobe  analysis  and  etch-pit  analysis.  The  MBE  growth 
conditions  employed  as  well  as  the  important  considerations  and 
procedures  pertaining  to  the  various  characterization  techniques 
mentioned  above  are  outlined  in  the  following  chapter. 

2.1  Molecular  Beam  Epitaxy 

Molecular  beam  epitaxy  is  in  essence  a vacuum  evaporation 
technique  but  with  two  important  differences -the  vacuum  is  ultrahigh 
vacuum  (UHV)  and  the  product  of  evaporation  is  a single  crystal.  In 
this  technique,  constituent  elements  and  dopants  are  directed  in  the 
form  of  thermal  molecular  or  atomic  beams  towards  a heated  single 
crystal  substrate  where  epitaxial  growth  occurs.  By  monitoring  the 
beam  fluxes  of  the  constituent  elements,  a precise  compositional 
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control  of  the  epitaxially  growing  layer  is  maintained.  High  purity 
source  elements  (in  our  case,  In,  Ga  and  As)  are  evaporated  from  high 
temperature  effusion  sources,  the  effusion  sources  being  arranged 
about  the  heated  crystalline  substrate  in  such  a way  as  to  ensure 
uniformity  both  of  composition  and  thickness.  Rotation  of  the 
substrate  during  growth  further  facilitates  the  attainment  of 
uniformity.  The  film  composition  is  varied  by  altering  the  source 
temperatures,  i.e.  the  temperature  of  the  furnaces,  thus  varying  the 
molecular  beam  fluxes.  Additional  control  over  the  growth  process  is 
achieved  by  virtue  of  mechanical  shutters  associated  with  each 
individual  source,  the  action  of  which  either  permits  or  prevents  the 
beam  from  reaching  the  substrate. 

With  regard  to  the  present  work,  the  following  characteristics  of 
MBE  make  it  an  ideal  technique  for  the  growth  of  multiple  layer 
InxGai-xAs/GaAs  structures. 

1.  In  general  the  growth  rate  is  low,  less  than  lpm/hr,  i.e.,  around  1 
monolayer/ sec,  which  allows  compositional  or  dopant  changes  to  be 
made  within  atomic  dimensions. 

2.  The  low  growth  temperatures  involved  minimize  interdiffusion 
effects. 

3.  Sequential  deposition  of  different  materials  is  possible  thus 
preserving  the  chemical  and  structural  integrity  of  the  various 
interfaces. 

4.  The  UHV  environment  allows  reflection  high  energy  electron 
diffraction  (RHEED)  analysis  to  be  employed  so  that  the  structural 
properties  of  the  epitaxial  layers  can  be  monitored  in  real-time. 
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In  this  study  a Varian  GEN  II  MBE  system  was  employed  for  the 
growth  of  InxGai-xAs  layers  on  GaAs  substrates,  a schematic  of  the 
growth  chamber  being  shown  in  Figure  2- 1 . The  GaAs  substrates  were 
typically  oriented  2°  off  the  (001)  plane  towards  the  (110)  plane  and 
prior  to  growth  the  substrates  were  exposed  to  a UV/ ozone 
environment  for  ~2  minutes.  This  procedure  resulted  in  the 
formation  of  a uniform  thin  oxide  film  on  the  substrate  surface  thus 
protecting  it  from  contaminants. 

After  being  exposed  to  UV/ozone,  the  substrates  were  attached 
to  molybdenum  blocks  using  indium,  and  the  molybdenum  blocks 
were  then  loaded  into  the  buffer  chamber  of  the  MBE  system  via  an 
entry/exit  chamber  which  were  pumped  down  to  pressures  of 
<5x1  O'®  torr  and  <lxl0"7  torr,  respectively.  The  molybdenum  block 
containing  the  substrate  was  then  transferred  into  the  growth 
chamber  where  the  oxide  was  desorbed  by  heating  the  substrate  at 
650°C  for  5 minutes  under  an  arsenic  overpressure. 

The  various  sources  used  in  the  growth  of  InxGai-xAs  layers 
included  one  elemental  In  source,  two  elemental  Ga  sources,  labeled 
Gal  and  Ga2,  and  one  tetramer  arsenic  (AS4)  source.  A 0.2(im  thick 
GaAs  buffer  layer  was  grown  homoepitaxially  on  the  substrate  prior  to 
the  growth  of  the  various  heterostructures  at  a growth  temperature  of 
600°  C. 

All  of  the  InxGa  i -xAs  layers  were  grown  employing  a beam 
equivalent  pressure  (BEP)  ratio,  i.e.  the  ratio  of  the  combined  fluxes  of 
the  group  III  elements  namely,  Ga  and  In,  to  the  flux  of  the  group  V 
element,  of  16-20.  The  fluxes  of  the  group  III  elemental  sources  and 
their  corresponding  source  temperatures  selected  to  provide  various 


28 


sources 


Figure  2-1.  Schematic  of  MBE  growth  chamber  configured  for  the 

growth  of  InxGai-xAs/GaAs  multiple  layer  systems.  The  Si 
and  Be  sources  are  n-  and  p-type  sources,  respectively. 
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InxGai-xAs  layer  compositions  are  provided  in  Table  2-1.  The 
substrate  temperatures,  i.e.  growth  temperatures,  were  optimized  [46] 
with  regard  to  growing  stoichiometric  InxGai-xAs  material  and  the 
various  temperatures  employed  are  given  in  Table  2- 1 . For  the  growth 
of  x<0.22  InxGai-xAs  layers  the  BEP  of  Ga  was  maintained  at 
5.33x1 0"7  while  the  BEP  of  In  was  altered  to  provide  the  requisite 
InxGai-xAs  layer  composition.  For  InxGai-xAs  layer  compositions  in 
the  range  0.22<x<0.4,  the  beam  equivalent  pressure  (BEP)  of  In  was 
maintained  at  2.26xl0-7  torr  while  the  Ga  pressure  was  altered  to 
provide  the  required  epilayer  composition.  Above  an  x value  of  0.4, 
the  BEP  of  In  was  maintained  at  a slightly  higher  value  of  3.33x1  O'7 
torr  in  order  to  provide  the  larger  amount  of  In  needed  for  these 
compositions  while  the  BEP  of  Ga  was  adjusted  to  provide  the 
required  composition. 

The  growth  of  multilayer  structures  consisting  of  layers  of 
InxGai-xAs  with  different  compositions  involves  growth  interrupts 
between  layers  consisting  of  two  different  compositions.  During  a 
growth  interrupt,  the  substrate  temperature  is  adjusted  to  correspond 
to  that  which  is  optimum  for  the  growth  of  the  subsequent  layer.  Also, 
the  In  source  is  shuttered  and  its  temperature  altered  to  provide  the 
requisite  beam  flux  for  the  growth  of  the  next  layer.  If  both  Ga  sources 
are  in  use,  then  they  are  both  shuttered  during  the  interrupt  and  their 
temperatures  are  adjusted  to  provide  the  required  beam  flux  for  the 
growth  of  the  succeeding  layer.  However,  if  only  one  of  the  Ga  sources 
is  active  during  the  growth  of  a layer  then  the  idle  or  shuttered  Ga 
source  temperature  is  altered  to  provide  the  optimum  layer 
composition  for  the  subsequent  layer  so  that  at  the  growth  interrupt. 
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Table  2-1.  Optimized  growth  conditions  employed  in  the  growth  of 
InxGai-xAs  layers  on  GaAs. 


Composition 

(x) 

Substrate 

temperature 

(°C) 

Group  III 
BEP 

(xl0'7torr) 

In  Ga 

BEP 
(V/ III) 
Ratio 

0.14 

520 

1.30 

5.33 

16-20 

0.22 

460 

2.26 

5.33 

16-20 

0.30 

420 

2.26 

3.52 

16-20 

0.35 

400 

2.26 

2.80 

16-20 

0.40 

380 

2.26 

2.26 

16-20 

0.45 

360 

3.33 

2.71 

16-20 

0.48 

360 

3.33 

2.40 

16-20 

0.50 

350 

3.33 

2.22 

16-20 

0.53 

340 

3.33 

1.97 

16-20 
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the  idle  Ga  source  is  activated  and  the  active  Ga  source  is  shuttered. 
This  cycle  is  repeated  at  each  growth  interrupt  encountered  during 
the  growth  of  a multilayer  structure.  Growth  of  the  various  layers  is 
monitored  using  RHEED  to  ensure  the  structural  integrity  of  the 
material.  While  a streaky  diffraction  pattern  designates  a material  of 
good  crystalline  quality  growing  in  a 2D  or  planar  growth  mode,  a 
spotty  pattern  represents  a rough  surface  caused  by  3D  or  island 
growth. 

After  growth,  the  epilayer/ substrate  structures  were  removed 
from  the  MBE  system,  detached  from  the  molybdenum  blocks  and 
extensively  analyzed  using  the  various  techniques  to  be  discussed 
below. 


2.2  Characterization  Techniques 

The  analytical  techniques  employed  in  this  study  included 
transmission  electron  microscopy,  etch-pit  analysis,  electron 
microprobe  analysis  and  high-resolution  x-ray  diffraction.  In  this 
section,  details  regarding  the  experimental  procedures  related  to  the 
various  characterization  tools  are  presented. 

2.2.1  Transmission  Electron  Microscopy 

The  various  InxGai-xAs/GaAs  heterostructures  grown  by  MBE,  as 
outlined  in  section  2.1,  were  analysed  by  cross-sectional  and  plan-view 
transmission  electron  microscopy  (XTEM  and  PTEM)  using  a JEOL 
200CX  transmission  electron  microscope.  XTEM  was  used  to  study 
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dislocation  evolution  in  the  heterostructures  and  PTEM  analysis  was 
used  to  determine  dislocation  densities  in  the  top  layers  of  stacked 
systems.  XTEM  and  PTEM  analyses  were  conducted  using  the  various 
reflections  available  on  (Oil)  and  (001)  zones,  respectively.  In  the 
following  subsections,  sample  preparation  techniques  for  XTEM  and 
PTEM  analyses  are  detailed. 

2. 2. 1.1  Cross-sectional  TEM  specimen  preparation  technique 

The  MBE  grown  samples  were  cut  into  10  mil  thick  strips  using 
a dicing  saw  and  the  top  surfaces,  i.e.  the  MBE  grown  epilayer 
surfaces,  of  two  of  the  strips  were  coated  with  M600  Bond  epoxy  resin 
and  brought  together  in  such  a way  that  these  coated  surfaces  were  in 
close  contact  with  each  other.  The  coated  strips  were  then  heated  for 
45  minutes  at  120°C  in  a conventional  oven  to  cure  and  thus  harden 
the  epoxy  resin.  Such  bonded  strips  were  polished  down  to  a 
thickness  of  about  50|im  using  a Struers  polishing  jig.  A copper  ring 
was  then  mounted  onto  the  thinned  specimen  to  provide  mechanical 
stability.  The  specimen  was  thinned  further  in  a Gatan  two-stage  ion 
mill  using  the  sputtering  action  of  two  Ar+  ion  guns  at  a gun  voltage  of 
5kV  and  a gun  current  of  0.5mA  until  a small  hole  was  obtained.  The 
regions  in  the  vicinity  of  the  hole  were  very  thin  and  thus  electron 
transparent  and  could  be  analysed  in  the  TEM. 
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2.2. 1.2  Plan-view  TEM  specimen  preparation  technique 

After  MBE  growth,  when  a sample  is  detached  from  the 
molybdenum  block,  residual  indium,  which  was  used  to  bond  the 
substrate  to  the  block  during  growth,  remains  on  the  substrate  side  of 
the  specimen.  Therefore,  the  first  step  towards  the  preparation  of  a 
PTEM  specimen  involves  polishing  the  substrate  side  of  the  sample 
until  the  indium  is  removed.  The  sample  is  then  cleaved  into  square 
pieces  in  such  a way  that  the  diagonal  length  of  each  square  is  about 
3mm.  A coat  of  wax  is  then  applied  to  the  top  surface  of  the  sample  to 
protect  it  during  the  etch.  The  sample  is  mounted  on  a holder  and 
etched  from  the  back  surface  in  a jet  etcher  using  a 1%  bromine- 
methanol  solution  until  an  electron  transparent  region  is  obtained. 

The  wax  coat  on  the  top  surface  is  then  removed  using  heptane  or 
trichloroethylene  and  the  specimen  is  ready  for  analysis. 

2.2.2  Etch  Pit  Technique 

Etch  pit  analysis  was  conducted  on  1pm  thick  InxGai-xAs  top 
layers  to  determine  the  threading  dislocation  density.  The  etchant 
solution  consisted  of  1ml  HF,  2ml  H2O,  8mg  AgNC>3  and  lg  003. 

Such  a solution  was  used  since  it  is  sensitive  to  dislocation  strain  fields 
and  thus  preferentially  etches  those  areas  leaving  etch  pits  behind.  A 
count  of  the  number  of  etch  pits  per  unit  area  provides  the  dislocation 
density.  By  measuring  the  thickness  of  the  epilayer  remaining  after  2, 
4,  8,  15  and  30  second-etches,  the  etch-rate  of  the  etching  solution 
was  determined  to  be  about  500A/second.  The  samples  were  etched 
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for  5 seconds  and  the  dislocation  density  was  determined  by  counting 
the  number  of  etch  pits  in  a lcmxO.lcm  area  using  a scanning 
electron  microscope. 

2.2.3  Electron  Microprobe  Analysis 

Electron  microprobe  analysis  was  used  to  determine  the 
composition,  x,  of  single  InxGai-xAs  layers  grown  on  GaAs  substrates. 
The  specimens  used  for  microprobe  analysis  consisted  of  InxGai-xAs 
layers  having  thicknesses  >l(im.  The  electron  microprobe  technique 
involves  the  use  of  an  electron  beam  to  probe  the  surface  of  the 
sample.  The  electrons  that  impinge  on  the  specimen  generate  x-rays 
which  are  characteristic  to  the  elements  comprising  the  material.  An 
electron  gun  voltage  of  only  6kV  was  used  in  order  to  limit  the  beam 
penetration  to  the  epilayer.  Integrated  x-ray  peak  intensities  of  the 
various  elements  were  evaluated  and  normalized  to  those  from 
reference  specimens,  i.e.  InAs  and  GaAs,  to  provide  the  composition  of 
the  InxGai-xAs  layer.  An  average  of  ten  measurements  across  the 
surface  of  the  specimen  were  made  to  evaluate  the  x value  of  the 
InxGai-xAs  epilayer.  The  error  involved  using  this  technique  was 
estimated  to  be  Ax<0.005. 

2.2.4  High  Resolution  X-Rav  Diffraction 

High  resolution  x-ray  diffraction  (HRXRD)  analysis  was 
performed  on  InxGai-xAs/GaAs  heterostructures  using  a Phillips  5- 
crystal  x-ray  diffractometer.  Since  HRXRD  is  a relatively  new  research 
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tool,  a brief  discussion  of  the  diffractometer  set-up  is  detailed  below. 

In  the  second  subsection  below,  the  optimization  procedures  required 
to  obtain  representative  HRXRD  rocking  curves  are  outlined  while  the 
various  steps  involved  in  the  measurement  of  strains  in  lattice- 
mismatched  systems,  such  as  the  InxGai-xAs/GaAs  material  system, 
are  discussed  in  the  third  subsection. 

2.2.4. 1 Phillips  high  resolution  x-rav  diffractometer 

The  Phillips  HRXRD  system  includes  an  x-ray  source,  a 
monochromator/  collimator,  a goniometer  and  a detector.  The  x-ray 
source  is  a Cu  target  anode  which  is  bombarded  by  accelerated 
electrons  created  at  a cathode  using  a generator  operating  at  40kV. 
Upon  such  bombardment,  the  Cu  target  emits  CuKa  x-ray  radiation  at  a 

characteristic  wavelength  of  1.54A..  The  x-ray  beam  thus  created  has  a 
broad  angular  and  wavelength  range  and  needs  to  be  monochromatized 
and  collimated  prior  to  its  impingement  upon  the  specimen.  The 
monochromator/ collimator  used  in  the  Phillips  system  consists  of  four 
highly  perfect  germanium  crystals  (Figure  2-2).  The  crystals  are  cut  in 
such  a way  that  either  their  (440)  or  (220)  planes  could  be  used  for 
diffraction.  The  crystals  are  aligned  so  as  to  reduce  spectral 
dispersion  and  peak  broadening  of  the  incident  x-ray  beam.  The  first 
crystal  (Figure  2-2)  is  set  to  the  correct  Bragg  angle  of  the  required 
diffracting  plane  (440  or  220)  with  respect  to  the  incident  Cukcc  x-ray 

beam.  The  schematic  in  Figure  2-2  shows  three  beam  paths 
corresponding  to  different  wavelength  components  of  the  dispersed 


36 


Figure  2-2.  A schematic  showing  the  alignment  of  the  four  Ge  crystals 
in  the  monochromator/collimator  section  of  the  Phillips 
HRXRD  system. 
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x-ray  beam.  All  three  components  are  diffracted  simultaneously  by  the 
first  crystal,  but  from  different  points  on  the  crystal  surface. 

The  second  crystal  is  fixed  in  position  such  that  it  satisfies  the 
Bragg  condition  for  the  median  (solid  line)  wavelength  component,  but 
as  shown  in  Figure  2-2  it  also  satisfies  the  Bragg  conditions  for  other 
wavelength  components.  Therefore,  there  is  no  reduction  in  the 
spectral  dispersion  of  the  x-ray  beam.  However,  the  reflectivity 
(intensity  of  the  diffraction  beam  about  the  Bragg  peak)  curve  of  the 
first  and  second  crystals  have  an  intrinsic  angular  width.  The 
superposition  of  the  reflectivity  curves  of  the  two  crystals  results  in  a 
reduction  in  the  angular  spread,  thus  collimating  the  x-ray  beam.  The 
third  crystal  is  set  in  a non-parallel  way,  and  only  a small  range  of 
wavelengths  can  satisfy  the  Bragg  condition  for  any  fixed  angle.  The 
range  of  wavelengths  passed  through  is  related  to  the  intrinsic  rocking 
curve  widths  of  the  reflections  used  at  the  third  crystal.  Due  to  this 
drastic  reduction  in  wavelength  range,  the  x-ray  beam  that  evolves  is 
monochromatic  in  nature.  The  fourth  and  third  crystals  are  aligned  in 
the  same  manner  as  crystals  one  and  two  and  hence  a combination  of 
crystals  3 and  4 serves  to  further  collimate  the  x-ray  beam.  The 
alignment  of  the  fourth  crystal  also  serves  in  reflecting  the  x-ray  beam 
back  onto  its  original  direction.  An  additional  advantage  of  such  a four 
crystal  monochromator/ collimator  set-up  is  the  ability  to  obtain 
rocking  curves  from  any  reflection  of  any  single  crystal  specimen. 

For  Ge  (440)  planes,  the  intrinsic  width  of  the  diffracted  beam  is 
about  5 arc  seconds  allowing  some  5%  of  the  characteristic  Ka  line  to 
pass.  Therefore,  the  use  of  (440)  reflections  in  the  four  Ge  crystals 
provides  an  extremely  monochromatic  beam,  however  with  very  low 
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intensity  (about  50,000  counts/second)  requiring  long  x-ray  scan 
times.  For  the  Ge  (220),  the  intrinsic  width  is  larger  than  the  (440) 
planes  and  is  about  12  arc  seconds,  thus  reducing  the  resolution  of  the 
diffractometer.  However,  the  intensity  of  the  exit  x-ray  beam  from 
(220)  planes  is  about  500,000  counts  per  second,  hence  requiring 
shorter  scan  times  are  required  for  the  aquisition  of  x-ray  data. 

The  x-ray  beam,  after  passing  through  the  monochromator/ 
collimator  section  impinges  on  a specimen  mounted  onto  the  sample 
stage  in  the  goniometer.  A schematic  of  the  goniometer  is  shown  in 
Figure  2-3.  The  goniometer  has  x,y  and  z cartesian  coordinate 
controls  for  the  adjustment  of  the  sample  stage  position.  The 
specimen  is  placed  on  the  sample  stage  and  the  z axis  of  the  stage  is 
adjusted  manually  in  such  a way  that  the  incident  x-ray  beam  and  the 
theta  axis  of  the  sample  lie  on  a plane  perpendicular  to  the  plane 
containing  the  sample  surface  normal  and  the  incident  beam.  The  x 
and  y axes  of  the  sample  stage  are  then  adjusted  via  computer  control 
so  that  the  x-ray  beam  is  incident  on  the  specimen-area  of  interest. 

The  goniometer  sample  stage  has  three  independent  degrees  of 
freedom  for  motion  in  spherical  coordinate  space.  The  sample  can  a) 
be  rotated  about  its  surface  normal  (z  axis)  over  a - 1 70°<R<+ 1 70° 
range,  b)  be  tilted  about  the  x axis  over  a -5°<Phi<+5°  range  and  c)  be 
rocked  about  the  y axis  over  a O°<20<18O°  range.  Furthermore,  the 
detector  can  move  independent  of  the  goniometer  over  a 340°  angular 
range.  The  angle  between  the  detector  and  the  sample  surface  is 
defined  as  20-to  where  co  is  the  angle  between  the  fixed  x-ray  source 
and  the  sample  surface.  One's  ability  to  rotate  and  tilt  the  specimen  is 
utilized  in  terms  of  the  optimization  procedures  to  be  discussed  in  the 
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Figure  2.3:  A schematic  illustrating  the  various  degrees  of  motion 
available  on  the  goniometer  sample  stage  used  for  x-ray 
rocking  curve  analysis  of  InxGai-xAs/GaAs  structures. 
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following  subsection  while  one's  ability  to  rock  the  specimen  and  move 
the  detector  independent  of  the  specimen  is  utilized  in  the  alignment 
of  the  diffractometer  for  the  procurement  of  x-ray  rocking  curve 
scans. 

As  illustrated  in  Figure  2-4,  two  types  of  HRXRD  rocking  curve 
scans  based  on  the  variation  of  co  and/or  20  were  used  in  this  work  as 
follows: 

a)  The  omega  scan:  where  the  detector  is  maintained  stationary  while 
the  specimen  rocks  over  a specified  angular  range  as  shown  in  Figure 
2-4a.  In  this  scenario,  20,  which  is  the  angle  between  the  incident 
beam  and  the  projected  diffracted  beam  that  reaches  the  detector 
remains  a constant  while  the  angle  between  the  specimen  surface  and 
the  x-ray  source,  co,  varies  over  a specified  range.  The  to  scan  is  useful 
for  the  determination  of  the  misorientation  of  a specimen  whose  Bragg 
condition  is  known.  The  20  value  is  fixed  to  satisfy  Bragg's  condition 
and  thus  at  a certain  value  of  to  (which  need  not  necessarily  equal  0)  an 
x-ray  peak  is  detected.  This  peak  position  can  be  used  to  determine 
the  specimen  geometry  which  will  be  discussed  in  detail  in  the 
following  section. 

b)  The  omega/2theta  scan:  where  the  range  of  to  is  specified. 

However,  the  specimen  is  rocked  over  the  specified  angular  range 
along  with  the  detector  which  sweeps  over  the  same  range.  The 
result  is  that  the  angular  variation  of  20  is  twice  that  of  co  (Figure  2-4b). 
Such  a scan  where  to  and  20  can  be  individually  specified  is  quite 
different  from  a 0/20  scan  performed  using  a conventional 
diffractometer  in  which  case  only  the  0 value  is  specified  and  which 
cannot  be  used  in  cases  where  to*0.  The  co/20  scan  is  used  when  the 
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Figure  2-4. 


Surface  normal 


(a) 


(b) 


Schematics  of  the  two  types  of  HRXRD  scans  employed  in 
the  present  study. 

a)  0)  scan  and  c)  co/26  scan. 
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specimen  is  composed  of  more  than  one  material  and  the  Bragg 
conditions  of  only  one  of  the  materials  is  known.  The  use  of  such  a 
scan  will  be  discussed  in  greater  detail  in  a later  section. 

2. 2. 4. 2 Optimization  procedures 

The  GaAs  substrates  used  for  growth  were  typically  miscut  at  2° 
off  the  (001)  plane  towards  the  (110)  plane.  Therefore  during  an  to 
scan,  the  position  of  the  GaAs  x-ray  peak  often  does  not  correspond  to 
its  real  value,  this  value  being  co=0  for  the  symmetric  (004)  peak.  A 
typical  substrate  miscut  geometry  is  shown  in  Figure  2-5.  As  seen  in 
the  figure,  as  the  specimen  is  rotated  about  its  [001]  normal  direction, 
the  position  of  the  peak  varies  in  the  range  to-Aa<to<co+Aa.  This 
variation  is  sinusoidal,  as  shown  in  Figure  2-6,  with  a periodicity  of 
180°,  the  difference  between  maximum  and  minimum  to  values 
corresponding  to  twice  the  misorientation. 

A layer  grown  on  a miscut  substrate  is  also  misoriented  as  seen 
in  Figure  2-5,  i.e.  the  surface  normal  does  not  correspond  to  the  [001] 
direction.  The  variation  of  peak  position,  co,  for  the  epilayer  is 
sinusoidal  just  as  it  is  for  the  substrate.  The  misorientation  of  the 
epilayer,  however,  is  not  identical  to  that  of  the  substrate  as  suggested 
by  the  differing  values  of  maxima  and  minima  for  epilayer  and 
substrate  indicated  in  Figure  2-6.  Such  mistilts  between  epilayer  and 
substrate  have  been  reported  for  other  epitaxial  systems  such  as 
InP/GaAs  [47,48],  AlxGai-xAs/  GaAs  [49]  and  ZnSe/GaAs[50-52].  The 
mistilts  have  been  explained  by  a simple  geometrical  model  assuming 
two-dimensional  growth  of  the  epilayer  which  starts  at  the  substrate 
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surface  normal 


Figure  2-5.  A schematic  showing  the  substrate  miscut,  Aa.  and  the 
epilayer  mistilt.  Ay,  with  respect  to  the  surface  normal. 
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Rotation  (R)  in  degrees 


Figure  2-6.  Plots  of  to  for  substrate  and  epilayer  as  a function  of 
specimen  rotation  angle.  R.  in  an  In0.lGa0.9As/GaAs 
structure  showing  a sinusoidal  variation  due  to  the 
substrate  miscut  (=2°)  and  an  epilayer  mistilt  (=0.5°). 
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steps  and  grows  in  such  a way  that  the  epilayer  surface  normal  (the 
growth  direction)  is  tilted  off  the  substrate  normal  and  towards  the 
(001)  direction.  This  explanation  appears  valid  for  the  InxGai-xAs/ 
GaAs  system  since,  as  shown  in  Figure  2-6,  the  misorientation  of  the 
epilayer  with  respect  to  the  [001]  direction  (deviation  from  optimal  co 
value)  is  less  than  that  between  the  substrate  and  the  [001]  direction. 

In  order  to  obtain  the  true  to  values  for  epilayer  and  substrate, 
the  rotation  angle  has  to  be  optimized.  Since  the  Bragg  angle  of  the 
epilayer  is  not  known,  to  scans  are  initially  conducted  at  a fixed 
substrate  Bragg  angle  (20  fixed)  and  a optimal  rotation  angle  for  the 
substrate  peak  is  obtained.  An  to/20  scan  is  then  performed  about  the 
substrate  peak  to  determine  the  Bragg  angle  of  the  epilayer.  This 
Bragg  condition  is  used  for  further  scans  to  optimize  the  epilayer  peak 
value.  In  order  to  obtain  accurate  values  for  to,  the  rotation  needs  to  be 
optimized  to  an  accuracy  of  1°. 

A second  method  for  the  determination  of  the  real  to  values 
involves  the  use  of  the  known  sinusoidal  nature  of  the  versus  R curves. 
The  sinusoidal  variation  implies  that  the  real  value  for  to  is  the  average 
of  any  two  to  values  separated  by  a rotation  angle  of  180°.  Therefore, 
rocking  curves  performed  at  two  rotation  angles  separated  by  180°, 
say  -90°  and  +90°,  would  be  required  to  determine  the  real  to  values 
for  epilayer  and  substrate. 

The  intensities  and  widths  of  the  various  x-ray  peaks  are 
critically  dependent  on  the  angular  orientation  of  the  specimen  with 
respect  to  the  x-y  plane,  i.e.  the  tilt  of  the  specimen  shown  in  Figure 
2-3.  Since  the  positions  of  the  peaks  are  not  affected  during 
optimization  of  the  tilt  angle,  co/20  scans  can  be  used  for  this  purpose. 
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The  advantage  of  using  such  scans  is  that  the  Bragg  conditions  for  both 
the  epilayer  and  substrate  peaks  are  satisfied  in  one  scan.  Therefore,  a 
number  of  x-ray  rocking  curve  experiments  are  performed  at  various 
tilt  angles  and  the  angle  that  provides  the  maximum  peak  intensity 
and  minimum  peak  width  is  chosen  as  the  optimum  value. 

The  terminology  for  the  peak  positions  used  in  the  subsequent 
section  represents  the  real  peak  positions  after  accounting  for  the 
miscut  of  the  substrate  and  the  mistilt  of  the  epilayer. 

2. 2. 4. 3 Strain  analysis  using  HRXRD 

The  peak  position,  to,  of  the  substrate  peak  at  a symmetric  plane, 
such  as  the  (004)  plane,  (symmetric  because  the  angle  of  incidence 
equals  the  angle  of  reflection  as  the  sample  surface  is  oriented  to  the 
(001)  plane)  can  be  used  to  determine  the  spacing,  d,  between  the 
(004)  planes  from  Bragg's  law  which  is  given  by 

nX  = 2dSin0 

where  X is  the  wavelength  of  the  incident  radiation  and  6 is  the  Bragg 

given  by  co=0  for  a symmetric  reflection.  The  angular  separation 
between  substrate  and  epilayer  peaks  for  the  symmetric  (004) 
reflection  can  be  used  to  determine  the  perpendicular  lattice- 
mismatch  between  epilayer  and  substrate  given  by 

^pi  — ^sub  N 

1 


v ^ub  y 


(®epi  ®sub)C0f®sub 
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where  a is  the  lattice  parameter  and  the  subscripts  "epi"  and  "sub" 
stand  for  the  epilayer  and  substrate,  respectively. 

In  general  one  must  distinguish  between  the  mismatch 
perpendicular  and  parallel  to  the  epilayer/ substrate  interface.  These 
mismatch  components  can  be  determined  by  using  HRXRD  rocking 
curves  from  asymmetric  lattice  planes  making  an  angle  <f>  with  the 

crystal  surface.  As  shown  in  Figure  2-7,  the  Bragg  condition  for  such  a 
plane  is  satisfied  at  two  different  co  angles  given  by 

wi  = 0 + <j> 

and 


C02  = 0 - <f> 


From  x-ray  rocking  curve  scans,  the  values  of  ©i  and  co2  can  be 
determined  for  both  epilayer  and  substrate.  By  solving  the 
aforementioned  simultaneous  equations,  the  values  for  6 and  <}>  for  both 
epilayer  and  substrate  can  be  determined.  A difference  in  lattice  plane 
spacing  between  epilayer  and  substrate  produces  a difference,  AG,  in 
the  Bragg  angle  for  a given  {hkl}  reflection,  whereas  a difference  in 
lattice  plane  orientation  corresponds  to  A<().  Assuming  tetragonal 
distortion  of  the  epilayer  lattice,  the  perpendicular  and  parallel 
lattice -mismatches  using  an  asymmetric  {hkl}  reflection  can  be 
evaluated  as  follows  [53,54]: 

Perpendicular  lattice-mismatch: 

^epl  — ^sub 

1 


V asub  , 


= A4>tan<J)sub  - A0cot0sub 
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detector 


(a) 


source 


(b) 


Figure  2-7.  Schematics  showing  two  conditions  at  which  Bragg's  law 

is  satisfied  for  asymmetric  planes  that  are  at  an  angle,  <J>,  to 
the  symmetric  (001)  plane. 

a)  co  1 = 0 - <)>  and  b)  co2  = 0 + <J) 
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Parallel  lattice -mismatch: 

= — A<t>cot<J>su5  - A0cot9sub 


^ ^pi  ^sub  ^ 


asub 


where  A0  and  A<>  are  given  by 


A0  = 0epi  - 6sub 


A<t>  = <t>epi  - <t>sub 


The  difference  between  the  perpendicular  and  parallel  lattice  - 
mismatches  is  related  to  the  deformed  epilayer  unit  cell.  Another 
method  to  determine  A0  and  A<|>  is  by  solving  simultaneously  the  peak 
separation  equations  from  rocking  curves  obtained  at  coi  and  a>2 
incidence,  respectively,  given  by 

Acoi  = A0  + A<(> 

Aco2  = A0  - A<> 

where  Acoi  and  Ao)2  are  the  peak  separations  at  the  two  angles  of 
incidence  of  the  x-ray  beam. 


CHAPTER  3 

STRAIN  RELIEF  IN  InxGai-xAs/GaAs  SINGLE  HETEROSTRUCTURES 


3.1  Introduction 

As  discussed  in  the  literature  review  section,  theoretical  and 
experimental  studies  conducted  on  the  InxGai-xAs/  GaAs 
heteroepitaxial  system  have  suggested  that  strain  associated  with  the 
lattice-mismatch  is  relieved  when  the  InxGai-xAs  epilayer  thickness 
exceeds  a critical  value  via  dislocation  formation  and  that  the  critical 
thickness  is  a function  of  the  InAs  mole  fraction  (x).  However,  as 
discussed  in  Section  1.2.3,  currently  published  theories  for  dislocation 
evolution  in  lattice-mismatched  systems  are  not  adequate  to  explain 
the  specifics  of  dislocation  evolution  in  the  InxGai-xAs/GaAs  material 
system.  In  this  chapter,  the  results  of  an  extensive  strain  relief  study 
are  detailed  and  a theory  is  presented  to  explain  dislocation  evolution 
in  the  In^ai-xAs/GaAs  heteroepitaxial  system. 

A set  of  InxGai-xAs  epilayers  was  grown  on  GaAs  with  a 0.5pm 
thick  GaAs  buffer  being  grown  prior  to  InxGai-xAs  growth  in  each 
case.  The  InxGai-xAs  epilayers  had  InAs  mole  fractions  of  x=0.07, 
0.14,  0.16,  0.18,  0.21,  0.28,  0.4  and  0.5  and  were  grown  to  a fixed 
thickness  of  1pm  which  is  well  above  the  critical  thickness  for  all  of 
the  compositions  studied.  Cross-sectional  transmission  electron 
microscopy  (XTEM)  analysis  was  performed  on  the  InxGai-xAs/GaAs 
single-layer  structures  in  order  to  observe  dislocation  evolution  as  a 


50 


51 


function  of  InxGai-xAs  composition.  The  XTEM  studies  were 
conducted  using  the  various  reflections  available  on  the  B[011]  zone 
axis. 

The  following  chapter  is  divided  into  three  sections.  The  first 
section  (3.2)  deals  with  the  experimental  results  from  XTEM  analysis 
of  the  various  heterostructures  while  the  second  section  (3.3) 
discusses  the  results  and  compares  them  to  published  reports. 

Finally,  a novel  theory  is  proposed  in  a third  section  (3.4)  based  on  a 
balance  of  forces  model  which  includes  a surface  image  force  term  to 
explain  dislocation  evolution  in  the  InxGai-xAs/GaAs  system. 

3.2  XTEM  analysis  of  InvGai  -vAs/GaAs  Single  Laver  Structures 

In  order  to  facilitate  an  explanation  of  the  observations  regarding 
dislocation  evolution,  the  InxGai-xAs/GaAs  structures  are  divided  into 
three  categories  as  follows: 

x < 0.18  Small  x structures 

0.18  < x < 0.28  Medium  x structures 
x > 0.28  Large  x structures 

Cross-sectional  TEM  analysis  of  the  small  x structures  indicated 
the  absence  essentially  of  threading  dislocations  in  the  InxGai-xAs 
epilayer  material  as  shown  in  Figure  3- la.  Threading  dislocations  are 
defined  in  this  work  as  dislocations  that  traverse  between  the 
luxGa  i -xAs / GaAs  interface  and  the  free  surface.  XTEM  analysis  of 
small  x In^a  i -xAs/  GaAs  heterostructures  revealed  a dark  band  at  the 
heterointerface  along  with  dislocations  in  the  GaAs  substrate  as  shown 
in  Figure  3- lb.  A tilting  experiment  on  such  a structure  showed  the 
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Figure  3-1.  Cross-sectional  TEM  micrographs  of  an  Ino.  l4Gao  86As/ 
GaAs  heterostructure. 

a)  A micrograph  showing  the  absence  of  dislocations  in 
the  epilayer;  b)  Dislocations  in  the  GaAs  substrate. 
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dark  band  to  be  a result  of  the  strain  field  from  misfit  dislocations  at 
the  heterointerface.  The  discrete  misfit  dislocations  can  be  clearly 
seen  in  the  XTEM  micrograph  of  Figure  3-2  obtained  from  a specimen 
that  was  severely  tilted  off  the  (Oil)  zone. 

Cross-sectional  TEM  analysis  of  medium  x structures  indicated 
the  presence  of  threading  dislocations  in  the  epilayer  as  well  as 
dislocations  in  the  GaAs  substrate  and  the  heterointerface  as  shown  in 
Figures  3-3a  and  3-3b.  Analysis  of  large  x (x>0.4)  InxGai-xAs/GaAs 
structures  showed  the  presence  of  threading  dislocations  in  the 
InxGai-xAs  epilayer  and  an  absence  of  dislocations  in  the  GaAs 
substrate  as  shown  in  Figures  3-4a  and  3-4b. 

Planar  TEM  analysis  was  performed  on  small  and  large  x 
structures  to  determine  threading  dislocation  densities  in  such 
epilayers  and  a comparison  of  the  dislocation  structure  in  a 
Ino.i4Gao.86As  (small  x)  layer  versus  that  in  a Ino.4Gao.6As  (large  x) 
layer  is  shown  in  Figure  3-5.  The  threading  dislocation  density  in  the 
Ino.4Gao.6As  epilayer  was  about  1010/cm2.  Since  the  PTEM 
specimen  of  the  Ino.i4Gao.86As  epilayer  showed  the  absence  of 
threading  dislocations,  an  estimate  of  the  threading  dislocation 
density  using  this  technique  was  not  possible.  The  etch-pit  technique 
was  therefore  required  to  determine  the  threading  dislocation  density 
in  such  small  x layers.  A scanning  electron  micrograph  of  an  etched 
specimen  is  shown  in  Figure  3-6.  The  calculated  threading  dislocation 
density  in  the  Ino.i4Gao.86As  epilayer  was  about  104/cm2. 
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Figure  3-2.  Transmission  electron  micrograph  of  a severely  tilted 
heterointerface  showing  discrete  misfit  dislocations. 
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Figure  3-3.  Cross-sectional  TEM  micrographs  of  a medium  x 

Ino.2lGao.79As  layer  showing  dislocation  formation  in 

both  the  epilayer  and  the  substrate. 

a)  Low  magnification  and  b)  High  magnification. 
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Figure  3.4:  Cross  sectional  TEM  micrograph  of  a large  x Ino.4Gao.6As 
layer. 

a)  Micrograph  showing  extensive  dislocation  formation  in 
the  epilayer;  b)  Micrograph  showing  a lack  of  dislocations 
in  the  substrate. 
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Figure  3.5:  Plan  view  TEM  micrographs  of  small  and  large  x layers. 

a)  Micrograph  of  an  Ino.  i4Gao.86As  layer  showing  a lack  of 
dislocations;  b)  Micrograph  of  a Ino.4Gao.6As  layer 
showing  an  extensive  network  of  dislocations. 
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3.3  Discussion  of  Results 

The  observation  of  dislocations  in  the  GaAs  substrate  in  small  x 
(x<0.18)  InxGai-xAs/GaAs  single-layer  structures  is  consistent  with  a 
previously  reported  study  of  Fitzgerald  et  al.  [28]  which  specifically 
concerned  Ino.l2Gao.88As/GaAs  samples.  The  force  acting  on 
dislocations  which  promotes  movement  into  the  GaAs  material  is  a 
repulsive  force  due  to  dislocations  having  the  same  Burgers  vector  at 
the  InxGai-xAs/GaAs  interface  while  the  major  opposing  force 
attracting  dislocations  back  to  the  heterointerface  is  an  interface 
image  force.  Thus  a balance  of  forces  would  cause  the  formation  of 
bowed  dislocations  or  dislocation  half-loops  in  the  GaAs  material. 

Such  bowed  dislocations  seemed  to  be  present  in  the  small  x 
structures  as  shown  in  Figures  3- la  and  3- lb.  However,  as  noted  in 
chapter  2,  the  absence  of  threading  dislocations  in  the  InxGai-xAs 
epilayers  (small  x)  was  not  well  explained  by  Fitzgerald  et  al.  [28]  who 
reasoned  that  the  larger  misfit  stress  in  the  epilayer  leads  to  the 
absence  of  threading  dislocations.  However,  such  a larger  misfit  stress 
in  the  epilayer  is  common  to  other  lattice-mismatched  systems  where 
an  absence  of  threading  dislocations  in  the  epilayer  is  not  observed. 

The  absence  of  threading  dislocations  in  the  epilayer  of  epilayer/ 
substrate  lattice-mismatched  systems  is  unique  to  the  InxGai-xAs/ 
GaAs  material  system  for  small  values  of  x.  This  unique  trait  can  be 
explained  on  the  basis  of  a comparison  of  the  high  temperature 
mechanical  properties  of  InxGai-xAs  versus  those  of  GaAs.  At  room 
temperature  InxGai-xAs  is  softer  (or  has  a lower  yield  strength)  than 
GaAs  [1].  However,  compression  tests  and  hardness  studies  conducted 
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Figure  3-6.  Scanning  electron  micrograph  of  an  etch-pit  in  a 

Ino.i4Gao.86As  epilayer.  Using  several  such  micrographs, 
the  dislocation  density  was  determined  to  be  lO^/cm^. 
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by  Bourret  et  al.  [55-57]  on  In  doped  GaAs  (x=0.001)  show  that  in 
comparison  to  GaAs,  the  yield  strength  of  the  In  doped  material 
decreases  more  gradually  with  increase  in  temperature.  The  results 
also  show  that  a critical  temperature  exists  at  which  the  yield 
strengths  of  In  doped  and  undoped  GaAs  material  are  identical,  the 
critical  temperature  being  found  to  be  390°C.  At  temperatures  greater 
than  the  critical  temperature,  the  In  doped  alloy  has  a greater  yield 
strength  than  that  of  the  undoped  GaAs  material  and  for  temperatures 
lower  than  the  critical  temperature,  undoped  GaAs  has  a greater  yield 
strength  than  In  doped  GaAs.  Also,  high  temperature  (780°C) 
compression  tests  conducted  by  Djemel  and  Castaing  [58]  on  In  doped 
GaAs  material  (x=0.004)  and  undoped  GaAs  show  that  the  yield 
strength  of  the  doped  GaAs  material  is  about  three  times  greater  than 
that  of  undoped  GaAs. 

Ehrenreich  and  Hirth  [59]  have  speculated  that  the  increase  in 
high  temperature  yield  strength  of  GaAs  upon  In  doping  is  due  to  solid 
solution  hardening.  In  the  case  of  metals,  single  solute  atoms  can 
produce  solution  hardening  [60].  For  GaAs,  by  contrast,  the 
corresponding  solute  entity  consists  of  5-atom  clusters  surrounded  by 
anions.  The  presence  of  a single  5-atom  InAs4  unit  in  GaAs  leads  to  a 
local  bond  extension  of  about  7%  (equal  to  the  lattice  parameter 
difference  between  InAs  and  GaAs)  or  equivalently  a volume  dilation  of 
21%.  Therefore,  this  unit  can  act  as  a strain  center  and  should  result 
in  substantial  solid  solution  strengthening.  The  assumption  that  the 
InAs4  unit  causes  all  of  the  observed  dilation  is  reasonable  since 
composition  dependence  studies  of  the  bond  lengths  on  InxGai-xAs 
alloys  [61]  show  that  the  As  sublattice,  while  retaining  a tetrahedral 
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configuration,  distorts  appreciably  in  such  a way  that  the  As  neighbors 
surrounding  the  Ga  atoms  contract  while  those  surrounding  In  atoms 
expand,  allowing  a volumetric  difference  of  21%  between  an  InAs4 
unit  and  its  surrounding  GaAs4  units. 

Guruswamy  et  al.  [621  conducted  high  temperature  (25-900°C) 
hardness  tests  on  In  doped  GaAs  alloys  (x=0.0025  and  0.01)  and  found 
that  the  high  temperature  hardness  of  the  Ino.OlGao.99As  alloy  is 
greater  than  that  of  the  Ino.0025Gao.9975As  alloy.  Since  the 
hardness  can  be  related  to  the  yield  stress,  it  was  determined  that  the 
yield  stress  of  the  Ino.OlGao.99As  alloy  is  greater  than  that  of  the 
In0.0025Gao.9975As  alloy  showing  an  increasing  trend  in  the  yield 
strength  with  increasing  x value.  Guruswamy  et  al.  [62J  noted  a soft- 
hard  transition  with  increasing  temperature  for  the  In  doped  alloys  in 
comparison  to  undoped  GaAs  but  could  not  obtain  an  exact  value  for 
the  transition  temperature  due  to  the  scatter  in  the  hardness  data. 
These  authors  speculated  that  the  transition  could  be  explained  based 
on  the  action  of  two  opposing  mechanisms,  a rate -controlling 
mechanism  related  to  the  intrinsic  lattice  resistance  or  Peierl's 
barrier  versus  a mechanism  for  dislocation  propagation.  At  low 
temperatures,  the  Peierl's  barrier  mechanism  is  operative  and  thus 
the  addition  of  In  could  soften  the  GaAs  material.  However,  as  the 
temperature  is  increased,  the  hardness  of  the  alloy  is  determined  by 
the  resistance  to  dislocation  glide  by  solute  atoms  or  InAs4  units  and 
thus  the  addition  of  In  increases  the  high  temperature  hardness  of 
GaAs. 

In  the  case  of  small  x InxGai_xAs  layers  grown  on  GaAs,  the 
smaller  yield  strength  of  GaAs  at  the  growth  temperature  allows  the 
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GaAs  substrate  to  yield  prior  to  the  InxGai-xAs  epilayer  upon 
application  of  strain  (such  as  lattice-mismatch  strain  in  the 
InxGai-xAs/GaAs  system)  and  hence  the  observation  of  dislocations  in 
the  GaAs  substrate  and  a lack  of  dislocations  in  the  Inj^Ga  i -xAs 
epilayer.  The  propagation  of  dislocations  into  the  GaAs  partially 
relieves  strain  in  addition  to  the  formation  of  misfit  dislocations  at  the 
heterointerface.  Also,  the  absence  of  dislocations  in  the  epilayer 
suggests  dislocation  nucleation  at  the  heterointerface. 

There  are  several  possible  theories  that  need  to  be  explored  in 
order  to  explain  the  evolution  of  dislocations  in  medium  and  large  x 
InjcGai-jcAs/GaAs  structures.  Based  on  the  argument  presented  above 
for  small  x InjcGai-jcAs/GaAs  structures,  i.e.  that  of  solid  solution 
hardening,  one  could  attempt  to  explain  the  presence  of  dislocations 
in  medium  and  large  x Inj^Gai^As  layers. 

It  is  possible,  for  instance,  that  the  solid  solution  hardening 
effect  due  to  In  addition  may  have  peaked  at  an  In^ai-jcAs 
composition  of  x=0.18.  This  peak  would  explain  the  presence  of 
threading  dislocations  in  the  InxGai-xAs  material  since  for 
InxGai-xAs/GaAs  structures  with  x>0.18,  the  yield  strength  of  the 
InjcGai-jcAs  epilayer  would  be  less  than  that  of  GaAs  allowing  the 
epilayer  to  yield  prior  to  the  GaAs  substrate.  However,  the  presence  of 
dislocations  in  the  GaAs  substrate  in  medium  x structures  suggests 
that  the  substrate  has  a comparable  yield  strength  to  that  of  the 
epilayer  at  the  corresponding  growth  temperature.  This  does  not 
seem  possible  since  the  presence  of  dislocations  in  both  the  substrate 
and  epilayer  occurs  over  a broad  range  of  compositions,  i.e. 
0.18<x<0.28.  In  actual  fact,  as  will  be  seen  in  chapter  5,  based  on  the 
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analysis  of  multilayer  structures,  the  yield  strength  of  InxGai-xAs 
appears  to  continually  increase  up  to  an  x value  of  ~0.5  at  elevated 
temperatures. 

A second  possible  explanation  for  the  presence  of  threading 
dislocations  in  medium  and  large  x InxGai-xAs  epilayers  could  involve 
a planar  (two-dimensional)  to  island  (three-dimensional)  growth  mode 
transition  at  InxGai-xAs  epilayer  compositions  >0.18.  In  this  case, 
during  growth,  the  InxGai-xAs  epilayer  would  be  deposited  in  the 
form  of  discrete  islands  which  would  grow  laterally  and  coalesce  giving 
rise  to  threading  dislocations  at  the  imperfect  island  junctions. 
Three-dimensional  growth  (island  growth)  can  be  observed,  in-situ,  by 
analytical  tools  such  as  reflection  high  energy  electron  diffraction 
(RHEED)  wherein  a streaky  pattern  would  signify  a smooth  surface  or 
planar  growth  and  a spotty  pattern  would  represent  a rough  surface  as 
in  the  case  of  island  growth.  However,  RHEED  analyses  on 
InxGai_xAs/GaAs  epilayers  performed  by  several  research  groups  [35- 
37]  have  shown  that  planar  growth  is  sustained  up  to  InxGai-xAs 
compositions  -0.5.  Also  the  growth  conditions  employed  in  this  work 
(as  mentioned  in  Chapter  2)  were  optimized  such  that  the  RHEED 
patterns  remained  streaky  during  growth.  Furthermore,  a 2D  to  3D 
growth  mode  transition  cannot  explain  the  presence  of  dislocations  in 
the  GaAs  substrate  in  the  medium  x InxGa  i -xAs/GaAs  structures. 
Therefore,  the  mechanism  of  2D-3D  growth  mode  transition  cannot  be 
utilized  to  explain  dislocation  evolution  in  medium  and  large  x 
InjcGai-jjAs/GaAs  structures. 

A third  possible  mechanism  by  which  threading  dislocation 
evolution  in  In^cGai^As  epilayers  could  be  explained  would  be  that  of 
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surface  nucleation  of  dislocations.  However,  as  indicated  in  the 
XTEM  micrographs  of  small  x structures,  nucleation  of  dislocations 
seems  to  occur  at  the  InxGai-xAs/GaAs  heterointerface,  suggesting 
that  interface  nucleation  is  energetically  favourable  when  compared  to 
surface  nucleation.  Therefore,  one  can  assume  that  the  first 
dislocations  to  form  in  medium  x InxGai-xAs  structures  would 
nucleate  at  the  heterointerface  when  the  critical  thickness  is  reached. 
These  dislocations  would  relieve  strain  in  the  epilayer  as  they  extend 
across  the  interface.  Consequently,  it  is  difficult  to  envision  surface 
nucleation  of  dislocations  given  the  larger  amount  of  energy  required 
to  nucleate  such  dislocations. 

In  the  following  section  a theory  based  on  a balance  of  forces 
acting  on  a dislocation  near  the  heterointerface  is  proposed  in  order 
to  explain  the  observed  dislocation  evolution  as  a function  of 
composition  in  the  InxGai-xAs/GaAs  material  system. 


3.4  A Force  Balance  Mechanism  for  Dislocation  Evolution 

As  mentioned  previously,  it  is  well  known  that  strain  relief 
occurs  by  the  formation  of  misfit  dislocations  above  a critical  thickness 
and  that  the  critical  thickness  decreases  with  increasing  lattice- 
mismatch  strain.  In  the  previous  section  it  was  shown  that  the 
presence  of  dislocations  in  the  substrate  material  as  opposed  to  the 
epilayer  can  only  be  explained  on  the  basis  of  dislocation  nucleation  at 
the  heterointerface.  Therefore,  at  the  onset  of  strain  relief,  i.e.  at  the 
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critical  thickness,  one  can  assume  that  misfit  dislocations  have  already 
nucleated  and  are  lying  in  the  plane  of  the  InxGai-xAs/GaAs  interface. 
The  force  balance  mechanism  to  be  discussed  below  will  only  be 
applicable  in  the  case  where  the  epilayer  is  harder  than  the  substrate 
at  the  growth  temperature  since,  if  the  converse  were  true,  surface 
nucleation  of  dislocations  would  become  a viable  option. 

In  general  terms,  several  forces  act  on  a misfit  dislocation 
segment  that  begins  to  propagate  into  an  epilayer  and,  depending  on 
the  net  force,  such  a segment  may  either  further  propagate  into  the 
epilayer  or  may  retract  towards  the  heterointerface.  Similarly, 
specific  forces  act  on  dislocations  that  may  propagate  into  the 
substrate. 

A dislocation  on  the  epilayer  side  of  the  heterointerface  would 
experience  a surface  image  force  due  to  its  proximity  to  the  free 
surface.  The  surface  image  force  (Fse)  would  attract  the  dislocation 
towards  the  free  surface  and  can  be  represented  as 

Geb2  (1  - vcos2p) 

Fse  47t(l  - v)d 


where  Ge  is  the  shear  modulus  of  the  epilayer  material,  d is  the 
thickness  of  the  epilayer,  b is  the  Burgers  vector,  v is  the  Poisson's 
ratio  and  P is  the  angle  between  the  dislocation  direction  and  the 
Burgers  vector.  A second  force  that  would  move  the  dislocation 
towards  the  surface  would  be  a repulsive  force  (Fde)  from  dislocations 
at  the  interface  that  have  the  same  Burgers  vector  and  is  given  by  [28] 
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where  the  first  term  represents  the  force  from  the  edge  component 
along  the  <11 1>  direction  (x  direction)  and  the  second  term 
represents  the  screw  component  resolved  along  x,  be  and  bs  are  the 
edge  and  screw  components  of  the  Burgers  vector  and  the  summation 
is  over  dislocations  at  the  interface. 

The  major  force  that  prevents  the  motion  of  dislocations 
towards  the  surface  is  the  force  from  the  elastic  strain  in  the  epilayer 
given  by 


where  e is  the  misfit  strain  in  the  epilayer  and  a is  the  lattice 
parameter  of  the  epilayer.  Other  forces  that  prevent  the  motion  of 
dislocations  towards  the  surface  include  the  Peierls  force  (Fpe).  an 
interface  image  force  (Fie)  and  a line  tension  force  (Fie)  which  are 
given  by  [63] 


where  0 is  the  angle  between  the  slip  plane  and  the  interface  plane 
and  R is  the  cutoff  radius  of  the  dislocation  which  is  typically  taken  as 
the  distance  to  the  interface  thus  making  the  expression  inaccurate 
when  near  the  interface  and 


where  d is  the  interplanar  spacing  between  (110)  planes, 
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Fie  = a 


GeP 
^ 4n  ( 1 - v)  > 


(1  - vcos2ft)  j'^xsinQ  j 


+ 2 


where  a is  a constant  based  on  the  geometiy  of  the  dislocation 
segment  and  p is  the  radius  of  curvature  of  the  segment. 

A dislocation  in  the  GaAs  substrate  would  experience  all  of  the 
aforementioned  forces  with  the  exception  of  the  surface  image  force 
[28].  Also,  all  of  the  forces  are  strongly  dependent  on  the  shear 
modulus  of  the  material  at  the  growth  temperature  [62]  and  hence, 
these  forces  will  be  weaker  in  the  GaAs  than  in  the  epilayer  material. 
Since  the  substrate  is  much  thicker  than  the  epilayer  material,  most  of 
the  elastic  strain  should  reside  in  the  epilayer  material.  Hence,  the 
force  from  the  elastic  strain  would  be  much  greater  in  the  epilayer 
than  in  the  substrate.  The  net  force  acting  on  a dislocation  in  the 
InxGai-xAs  epilayer  allowing  it  to  move  towards  the  free  surface  can 
be  expressed  as 


Fepi  - Fse  + Fde  - Fpe  - Fee  ■ Fie  - Fie 


while  the  net  force  acting  on  a dislocation  in  the  GaAs  allowing  it  to 
traverse  further  into  the  substrate  is  given  by 

Fsub  = Fds  - FpS  - Fes  - FiS  - Fis 


where  the  subscript  s stands  for  the  force  acting  on  a dislocation  in 
the  GaAs.  Of  the  various  forces  acting  on  a dislocation,  the  surface 
image  force  and  the  elastic  strain  force  are  the  ones  that  are 
extremely  dependent  on  InxGai-xAs  composition.  As  the  composition 
of  InxGai-xAs,  i.e.  the  x value,  increases,  the  critical  thickness 
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decreases  due  to  the  increasing  lattice-mismatch  strain.  Therefore, 
the  surface  image  force  increases  with  increasing  InxGai-xAs 
composition.  However,  the  major  force  that  prevents  a dislocation 
from  traversing  to  the  surface,  i.e.  the  elastic  strain  force  (Fee),  also 
increases  with  increasing  InxGai-xAs  composition.  A ratio  of  the 
surface  image  force  to  the  elastic  strain  force  as  a function  of  epilayer 
composition,  x,  is  plotted  in  Figure  3-7.  The  critical  thicknesses  used 
to  calculate  the  surface  image  force  values  were  obtained  from 
Matthews  equation  to  predict  the  critical  thickness  given  by 


b 

87te(  1 + v) 


1 + In 


where  e is  the  misfit  strain.  The  plot  in  Figure  3-7  indicates  that  for 
small  x structures,  the  surface  image  force  is  small  due  to  the  large 
values  of  critical  thickness.  Also,  since  Fee»Fes.  Fepi  is  much  smaller 
than  Fsub-  Therefore,  dislocations  in  the  small  x structures  should 
appear  only  in  the  GaAs  as,  in  fact,  has  been  observed.  At  some 
"critical  composition"  the  surface  image  force  would  become 
comparable  to  the  elastic  strain  force,  thus  allowing  dislocations  to 
traverse  to  the  free  surface.  This  composition,  as  narrowed  down  by 
XTEM  analysis,  appears  to  be  at  x=0.18. 

In  medium  x structures,  Fse  becomes  significantly  larger  and 
hence  Fepi  can  approximate  FSub-  Thus,  the  observation  of 
dislocations  in  both  the  epilayer  and  the  GaAs  in  this  case.  Finally,  for 
very  large  values  of  FSe,  as  in  the  case  of  large  x structures,  Fepi  will 
become  much  greater  than  Fsub  and  hence  all  of  the  dislocations  will 
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Figure  3-7.  A plot  of  the  ratio  of  the  surface  image  force  (Fse)  to  the 
misfit  force  (Fee)  versus  InxGai-xAs  composition  showing 
an  increase  in  the  surface  image  force  with  increasing 
composition. 
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be  able  to  travel  to  the  free  surface  leaving  the  GaAs  dislocation  free  as 
observed  in  this  work. 


3.5  Summary 

In  summary,  it  has  been  shown  that  a "critical  composition"  of 
x=0.18  exists  in  the  InxGai-xAs/GaAs  system  up  to  which  InxGai-xAs 
layers  can  be  grown  directly  on  GaAs  without  the  evolution  of 
threading  dislocations  in  the  epilayer.  This  trait  of  the  InxGai^As/ 
GaAs  material  system  has  been  explained  by  assuming  that  InxGa  i -xAs 
has  a greater  yield  strength  than  that  of  GaAs  at  the  growth 
temperature  which  causes  the  substrate  to  yield  preferentially  and 
hence  relax  the  misfit  strain.  Also,  dislocation  evolution  in  the 
InjcGai-jcAs/GaAs  system  as  a function  of  x has  been  explained  by 
assuming  that  a balance  of  forces  acting  on  misfit  dislocations  near  the 
heterointerface  is  attained.  The  magnitude  of  the  various  forces, 
which  are  very  much  alloy  composition  dependent,  determine 
whether  or  not  dislocations  propagate  in  the  epilayer  or  the  substrate, 
or  indeed,  both  the  epilayer  and  the  substrate. 


CHAPTER  4 

STRAIN  RELIEF  IN  InxGai-xAs-BASED  STRAINED  LAYER 
SUPERLATTICE  SYSTEMS 

4.1  Introduction 

In  the  previous  chapter  it  was  shown  that  low-dislocation 
density  (104/cm2)  InxGai-xAs/GaAs  epilayers  can  be  grown  directly 
on  GaAs  for  x values  <0.18.  However,  epilayers  having  x values  greater 
than  0.18  were  shown  to  contain  a high  density  of  threading 
dislocations.  The  work  to  be  described  in  the  following  chapter 
concerned  the  use  of  InjcGai-xAs-based  strained  layer  superlattices 
(SLSs)  which  were  investigated  with  a view  towards  reducing  the 
threading  dislocation  density  in  In^ai-xAs/GaAs  epilayers  having  x 
values  greater  than  0.18. 

4.2  Single  Strained  Laver  Superlattice  Buffer  Structures 

As  discussed  in  the  literature  review  section,  strained  layer 
superlattices  have  been  shown  to  be  somewhat  effective  in  blocking 
dislocations  from  propagating  into  epilayers  in  lattice-mismatched 
systems  in  which  case  dislocation  bending  by  strain  fields  appears  to 
be  the  primary  factor  responsible  for  dislocation  blocking.. 
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In  order  to  determine  the  effectiveness  of  such  SLSs  in  terms  of 
reducing  the  dislocation  density,  for  example,  in  an  Ino.4Gao.6As/ 
GaAs  epilayer,  a SLS  consisting  of  100  periods  of  10A  thick  alternate 
layers  of  GaAs  and  Ino.8Gao.2As  was  grown  on  GaAs  and  a 1pm  thick 
Ino.4Gao.6As  layer  was  grown  on  the  SLS.  The  x=0.8  composition  was 
chosen  in  this  case  in  order  that  the  average  composition  of  the  SLS 
match  that  of  the  epilayer.  XTEM  analysis  of  the  structure  showed 
that  dislocations  which  formed  at  the  GaAs/SLS  heterointerface 
threaded  through  the  SLS  and  the  InxGai-xAs  epilayer  to  the  free 
surface  (Figure  4-1).  It  was  therefore  concluded  that  this  particular 
SLS  was  ineffective  as  a dislocation  blocker.  The  ineffectiveness  of  the 
SLS  could  be  due  to  the  large  lattice-mismatch  strain  of  3%  between 
substrate  and  SLS  and/or  the  large  mismatch  strain  of  6%  between  the 
individual  layers  in  the  SLS.  In  a second  experiment,  a single 
Ino.3Gao.7As/  Ino.5Gao.5As  SLS  was  grown  directly  on  GaAs  so  that 
the  mismatch  between  individual  layers  in  the  SLS  would  be 
significantly  reduced  when  compared  to  the  previously  studied  SLS, 
while  the  mismatch  between  the  substrate  and  SLS  would  be  the  same 
as  in  the  previous  case  (3%).  The  thickness  of  individual  layers  was 
maintained  below  the  critical  thickness  at  a value  of  100A  and  a 1pm 
thick  Ino.4Gao.6As  layer  was  grown  on  top  of  the  SLS.  XTEM  analysis 
of  the  heterostructure  showed  dislocations  propagating  through  the 
superlattice  to  the  top  Ino.4Gao.6As  layer  (Figure  4-2a).  A closer 
observation  of  the  SLS  (Figure  4-2b)  indicated  that  the  dislocations 
were  not  at  all  bent  by  the  stresses  within  the  individual  SLS  layers. 

The  ineffectiveness  of  such  SLSs  was  attributed  to  the  large  mismatch 
strain  existing  between  the  substrate  and  the  SLS  and,  in  order  to 
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Figure  4-1.  Cross-sectional  TEM  micrograph  showing  a single  SLS 
incorporated  between  the  GaAs  substrate  and  an 
Ino.4Gao.6As  layer.  The  SLS  is  not  resolved  due  to  the 
extremely  small  thicknesses  of  individual  layers.  Extensive 
threading  dislocation  formation  is  noticed  in  the  epilayer. 
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Figure  4-2.  Cross-sectional  TEM  micrographs  of  a single  SLS 
structure  designed  to  block  dislocations, 
a)  Micrograph  of  the  entire  structure;  b)  High- 
magnification  of  the  SLS  showing  its  ineffectiveness  in 
blocking  dislocations. 
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achieve  low  dislocation  density  large  x InxGai-xAs  layers,  it  was 
decided  that  perhaps  a multistage  buffer  system  would  be  necessary 
wherein  strain  relief  would  occur  on  a gradual  basis  upon  growth  of 
each  stage.  In  the  following  section  the  design  and  analysis  of  such  a 
strain  relief  structure  is  detailed. 

4.3  Design  and  Analysis  of  Multistage  Strain  Relief  Systems  Containing 

Strained  Laver  Superlattices 

The  buffer  system  in  this  work  comprised  stacked  stages  which 
were  designed  in  such  a way  as  to  promote  strain  relief  and  dislocation 
blocking  in  underlying  layers  of  the  stack  leaving  the  top  layer 
relatively  dislocation-free.  XTEM  analysis  was  conducted  after  each 
stage  of  the  stacked  structure  was  grown  and  design  of  each 
subsequent  stage  was  based  on  this  analysis.  The  additional  stage  was 
incorporated  in  the  multistage  heterostructure  which  was  then  grown 
anew. 

A 0.2|im  thick  Ino.i4Gao.86As  layer  was  chosen  for  the  first 
stage  of  the  multistage  system  since  such  a composition  lies  below  the 
critical  composition  of  x=0.18.  Strain  relief  at  this  stage  occurs  by  the 
evolution  of  dislocations  at  the  Ino.i4Gao.86As/GaAs  heterointerface 
and  in  the  GaAs  substrate  as  explained  in  the  previous  chapter. 

A 0.2|im  thick  Ino.27Gao.73As  layer  was  grown  at  the  second 
stage  so  that  the  difference  in  compositions  between  the  first  and 
second  stages  would  be  below  the  critical  composition.  However, 

XTEM  analysis  of  the  two-stage  system  revealed  a high  density  of 
threading  dislocations  in  the  top  layer  (Figure  4-3).  At  this  point  it 
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Figure  4-3.  Cross-sectional  TEM  micrograph  of  a two-layer  structure 
showing  extensive  threading  dislocation  generation  in  the 
top  Ino.27Gao.73As  layer. 
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was  decided  that  a SLS  would  be  necessary  between  the 
Ino.14Gao.86As  and  Ino.27Gao.73As  layers  to  block  dislocations  from 
propagating  into  the  top  layer.  The  compositions  of  the  alternating 
layers  in  the  SLS  were  chosen  to  be  x=0.14  and  x=0.4  such  that  the 
average  composition  of  the  SLS  would  equal  the  top  layer  composition. 
The  thicknesses  of  the  layers  were  chosen  to  be  20A  which  is  less 
than  the  critical  thickness  and  therefore  dislocation  nucleation  within 
the  SLS  structure  would  not  be  a concern.  An  XTEM  micrograph  of 
the  two  stage  structure  consisting  of  a single  InxGai-xAs  layer  at  stage 
1 and  a 50  period  SLS  at  stage  2 is  shown  in  Figure  4-4a.  As  can  be 
seen  from  the  figure,  dislocation  formation  in  the  Ino.i4Gao.86As 
layer  is  evident  along  with  the  presence  of  a few  dislocations  in  the 
SLS.  Dislocations  in  the  SLS  appeared  to  be  effectively  bent  by  the 
stresses  present  in  the  SLS  layers  as  evidenced  by  the  lack  of 
dislocations  in  the  layers  near  the  surface.  A 0.2pm  Ino.27Gao.73As 
layer  was  chosen  for  the  third  stage  since  it  should  be  lattice-matched 
to  the  underlying  SLS.  XTEM  analysis  of  the  three-stage  structure 
showed  a lack  of  threading  dislocations  in  the  top  layer  (Figure  4-4b). 
A comparison  of  this  three-stage  structure  to  the  structure  shown  in 
Figure  4-3  clearly  shows  the  advantage  of  using  a SLS  to  block 
dislocations.  HRXRD  analysis  was  performed  to  study  lattice-matching 
between  the  second  and  third  stages  of  the  structure  shown  in  Figure 
4-4.  Typical  rocking  curves  at  a (444)  reflection  obtained  from  two- 
and  three-stage  structures  (Figure  4-5)  indicated  two  epilayer  peaks 
for  both  structures.  In  the  two-stage  structure  these  two  epilayer 
peaks  correspond  to  stages  1 and  2 of  the  heterostructure, 
respectively.  The  peak  representing  the  SLS  at  stage  2 corresponds 
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Figure  4-4.  Cross-sectional  TEM  micrographs  of  two-  and  three-stage 
structures  showing  effective  dislocation  blocking  by  the 
SLS  at  stage  2 resulting  in  a low  dislocation  density 
Ino.27Gao.73As  layer. 

a)  Two-stage  structure;  b)  Three-stage  structure. 
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Figure  4-5.  X-ray  rocking  curves  obtained  from  the  co"  condition  of 
(444)  planes. 

a)  Rocking  curve  from  a two-stage  structure;  b)  Rocking 
curve  from  a three-stage  structure  showing  an  overlap  of 
peaks  from  stages  2 and  3. 
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to  its  Oth  order  peak.  Higher  order,  or  satellite,  peaks  are  not  seen  in 
the  rocking  curve  shown  in  Figure  4-5a  due  to  the  limited  angular 
scan  range  (Aco)  employed.  A HRXRD  analysis  of  the  three -stage 
structure  revealed  that  stages  2 and  3 are  represented  by  a single 
epilayer  peak  (Figure  4-5b).  Therefore,  one  can  conclude  that  the  two 
stages  are  lattice-matched  in  both  the  in-plane,  or  parallel,  and 
perpendicular  directions. 

A second  SLS  was  chosen  for  use  at  the  fourth  stage  of  the 
multistage  structure  and  designed  such  that  the  difference  between 
the  average  composition  of  the  SLS  and  that  of  the  underlying  layer 
(x=0.27)  was  below  the  critical  composition  of  x=0.18.  The 
compositions  of  the  individual  InxGai-xAs  layers  comprising  the  SLS 
were  chosen  to  be  x=0.27  and  x=0.53,  respectively,  giving  an  average 
SLS  composition  of  x=0.4.  The  thicknesses  of  the  individual  layers 
were  20A  each  and  there  were  50  periods  forming  the  SLS.  XTEM 
analysis  of  the  four-stage  structure  showed  dislocation  evolution  in  the 
underlying  third  stage  (Figure  4-6a).  Dislocations  were  also  observed 
in  the  SLS  grown  at  the  fourth  stage.  Most  of  the  dislocations  in  the 
SLS  seemed  to  be  blocked  and  diverted  toward  the  edges  of  the 
specimen.  However,  a few  dislocations  appeared  to  traverse  through 
the  SLS  to  the  free  surface.  These  dislocations  appeared  to  be 
nucleating  within  the  SLS  rather  than  at  the  Stage  3/  Stage  4 
interface.  Since  the  mismatch  between  alternating  layers  in  the  SLSs 
at  both  stages  2 and  4 are  equal  (lattice  parameter  difference,  Aa,  is 
constant  since  compositional  difference,  Ax,  is  constant  at  a value 
Ax=0.26),  the  nucleation  of  dislocations  or  yielding  within  the  SLS  at 
stage  4 could  be  explained  assuming  that  the  yield  strength  of  one  of 
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Figure  4-6.  Cross-sectional  TEM  micrographs  of  four-  and  five-stage 
structures. 

a)  Micrograph  of  a four-stage  structure  showing  strain 
relief  in  stages  1-3.  The  fourth  stage  SLS  shows  a few 
dislocations  that  appear  to  propagate  to  the  free  surface. 
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Figure  4-6— continued 

b)  Micrograph  of  a five-stage  structure  showing  a lack  of 
dislocations  in  the  top  Ino.4Gao.6As  layer. 
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the  layers  in  the  SLS  is  smaller  than  that  of  the  single  InxGai-xAs 
layer  at  stage  3.  The  yield  strength  as  a function  of  InxGai-xAs 
composition  will  be  discussed  in  greater  detail  in  Chapter  5. 

A 1pm  thick  single  InxGai-xAs  layer  with  a composition  of  x=0.4 
was  chosen  for  the  fifth  stage  due  to  its  match  to  the  average  lattice 
parameter  of  the  SLS  at  stage  4.  XTEM  analysis  of  the  structure 
revealed  that  the  top  InxGai-xAs  layer  was  relatively  dislocation-free 
(Figure  4-6b)  while  strain  relaxation  occurred  by  dislocation  evolution 
in  the  underlying  stages,  i.e.  stages  1-4.  HRXRD  analyses  of  four-  and 
five-stage  structures  showed  the  presence  of  four  epilayer  peaks  in  the 
respective  x-ray  rocking  curves  (Figure  4-7).  The  presence  of  four 
epilayer  peaks  corresponding  to  the  four  stages  in  the  four-stage 
structure  (Figure  4-7a)  suggests  that  the  superlattice  at  stage  2 
relaxed  to  a different  extent  when  compared  to  the  single  InxGai-xAs 
layer  at  stage  3 upon  growth  of  the  fourth  stage,  since,  prior  to  the 
growth  of  the  fourth  stage,  stages  2 and  3 were  lattice-matched 
(Figure  4-5b).  Such  a relaxation  difference  between  stages  2 and  3 can 
be  expected  since  the  SLS  is  made  up  of  two  different  layers  (in  terms 
of  composition)  whose  mechanical  properties  are  very  different. 
Therefore,  one  half  of  the  SLS  could  yield  at  a vastly  different  rate  in 
comparison  to  the  other  half  and  consequently,  the  net  strain 
relaxation  of  the  SLS  would  be  different  when  compared  to  that  of  the 
single  layer  stage  whose  composition  and  hence  properties  are 
different  from  those  of  the  layers  in  the  SLS. 

Plan-view  TEM  analysis  of  the  five-stage  structure  (Figure  4-8) 
revealed  that  the  threading  dislocation  density  in  the  Ino.4Gao.6As  top 
layer  was  about  10®/ cm2.  Such  a defect  density  has  proved  to  be 
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Figure  4-7.  X-ray  rocking  curves  obtained  using  the  to-  condition  on 
(444)  planes. 

a)  Rocking  curve  from  a four-stage  structure  showing 
discreet  peaks  representing  stages  2 and  3.  respectively; 

b)  Rocking  curve  from  a five-stage  structure  showing  an 
overlap  of  peaks  representing  stages  4 and  5. 
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Figure  4-8.  Plan-view  TEM  micrograph  of  the  Ino.4Gao.6As  top  layer 
showing  a few  threading  dislocations.  The  dislocation 
density  was  determined  as  being  10^/cm2. 
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sufficiently  low  for  the  production  of  good  quality  Ino.4Gao.6As/GaAs 
photodetectors  [64].  However,  a much  lower  dislocation  density 
(<105/cm2)  would  probably  be  required  for  the  production  of  more 
defect  sensitive  devices  such  as  LEDs  and  diode  lasers.  A close 
examination  of  the  stage  4 superlattice  in  the  five -stage  structure 
(Figure  4-9)  showed  that  some  dislocations  were  unaffected  by  the 
superlattice  and  were  able  to  traverse  to  the  top  layer  thus 
contributing  to  the  threading  dislocation  density  observed  by  PTEM. 

In  an  effort  to  design  more  effective  dislocation-blocking  SLSs,  the 
thicknesses  of  the  individual  layers  in  the  SLS  at  the  fourth  stage  was 
altered  to  10A,  50A  and  100A  in  three  separately  grown  five-stage 
structures  so  that  the  stresses  in  the  SLS  layers  would  be  different. 
The  compositions  of  the  individual  layers  comprising  the  SLS  and  the 
total  number  of  periods  were  maintained  the  same  as  before.  XTEM 
analysis  of  the  SLSs  indicated  no  change  with  regard  to  the 
effectiveness  of  the  various  superlattices  in  blocking  dislocations. 

In  order  to  study  the  compositional  range  over  which  such  low 
dislocation  density  InxGai-xAs  layers  could  be  obtained  by  the  use  of 
the  multistage  strain  relief  buffer  system,  a seven-stage  structure  was 
designed  wherein  the  sixth  stage  consisted  of  a SLS  and  the  seventh 
stage  consisted  of  a 1pm  thick  Ino.48Gao.52As  layer  that  was  lattice- 
matched  to  the  underlying  superlattice.  The  SLS  comprised  50 
periods  of  20A  thick  alternating  Ino.4Gao.6As  and  Ino.56Gao.44As 
layers.  XTEM  analysis  of  the  seven-stage  structure  showed  extensive 
threading  dislocation  evolution  in  the  top  Ino.48Gao.52As  layer  (Figure 
4-10).  PTEM  analysis  of  the  top  layer  revealed  that  the  threading 
dislocation  density  was  about  108/cm2  (Figure  4-11). 
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Figure  4-9.  Cross-sectional  TEM  micrograph  of  the  fourth  stage  SLS 
in  the  five-stage  structure  showing  a dislocation 
propagating  through  the  superlattice  to  the  top  layer 
without  being  affected  by  the  stresses  in  the  SLS. 
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Figure  4-10.  Cross-sectional  TEM  micrograph  of  a seven-stage 
structure  showing  extensive  threading  dislocation 
formation  in  the  top  Ino.48Gao.52As  layer. 
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Figure  4-11.  Plan-view  TEM  micrograph  of  a Ino.48Gao.52As  top  layer 
showing  a high  density  of  threading  dislocations  despite 
the  use  of  a strain  relief  buffer. 
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4.4  Summary  and  Conclusions 

A multistage  strain  relief  buffer  structure  involving  several  SLSs 
as  opposed  to  a single  SLS  was  shown  to  be  fairly  effective  in  terms  of 
reducing  the  density  of  dislocations  threading  in  InxGai-xAs  layers 
having  x values  >0.18.  The  smallest  dislocation  density  obtained  by 
use  of  a multistage  buffer  in  the  case  of  Ino.4Gao.6As  material  was 
~10®/cm2.  However,  attempts  at  further  dislocation  density 
reduction  by  varying  the  SLS  structure  proved  unsuccessful,  which  was 
probably  due  to  dislocation  generation  within  the  SLSs  themselves.  In 
terms  of  growing  relatively  low  dislocation  density  InxGai-xAs  on  GaAs 
employing  SLSs,  a maximum  In^cGa  1 -xAs  composition  of  x=0.4  was 
determined.  Clearly,  a different  method  of  accomodating  strain  relief 
would  be  necessary  to  provide  low  dislocation  density  large  x 
InxGai-xAs  epilayers  on  GaAs. 


CHAPTER  5 

STRAIN  RELIEF  IN  STACKED  InxGai-xAs  MULTILAYER  SYSTEMS 
NOT  INVOLVING  STRAINED  LAYER  SUPERLATTICES 

5.1  Introduction 

As  concluded  in  the  previous  chapter,  a strain  relief  system 
other  than  one  involving  SLSs  would  be  required  in  order  to  extend 
the  composition  range  over  which  low-dislocation  density  (<l()5/cm2) 
InxGai-xAs  could  be  grown.  With  such  an  objective  in  mind,  it  was 
decided  that  in  order  to  effectively  utilize  the  strain  relief  mechanism 
proposed  in  chapter  3 toward  the  development  of  multilayer 
InxGai-xAs  structures,  a quantitative  understanding  of  the  extent  of 
strain  relaxation  occurring  in  each  layer  of  a multilayer  system  would 
be  vital.  To  this  end,  the  extent  of  lattice  relaxation  occurring  in 
single  InxGai-xAs/GaAs  heterostructures  as  functions  of  composition 
and  layer  thickness  was  first  of  all  studied  and  the  results  of  this  work 
are  detailed  in  section  5.2  below.  This  work  was  then  extended  to 
stacked  multilayer  InxGai-xAs  systems,  the  extent  of  lattice  relaxation 
occurring  in  each  layer  of  a stacked  system  being  determined  as  a 
function  of  the  number  of  layers  forming  the  stacked  structures. 
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5.2  Extent  of  Residual  Strain  in  Single  In-rGai  -vAs/GaAs 

Heterostructures 

In  order  to  study  the  extent  of  strain  relaxation  in  single 
InxGai-xAs  layers  grown  on  GaAs,  a variety  of  layers  were  grown 
having  various  compositions  (as  determined  by  electron  microprobe 
analysis)  and  thicknesses  and  high  resolution  x-ray  diffraction 
(HRXRD)  analysis  was  performed  on  the  samples.  All  of  the  x-ray 
rocking  curves  were  obtained  after  optimizing  the  tilt  of  the  specimen 
for  maximum  substrate  peak  intensity  as  explained  in  Chapter  2. 

Since  the  GaAs  substrates  were  oriented  2°  off  (001)  towards  [1101, 
the  x-ray  peak  position  (co)  varied  sinusoidally  with  speciman  rotation 
angle  (R)  with  a periodicity  of  360°.  In  order  to  determine  the  real 
peak  positions,  rocking  curves  were  obtained  at  two  specimen  rotation 
angles  separated  by  180°  and  the  average  of  the  peak  positions  from 
the  two  scans  was  taken  as  the  true  value  of  to.  HRXRD  analysis  was 
conducted  on  asymmetric  (444)  planes  and  the  perpendicular  and 
parallel  lattice-mismatches  were  evaluated  with  respect  to  the  GaAs 
substrate  assuming  tetragonal  distortion. 

An  example  of  residual  strain  analysis  concerning  a 0.5pm  thick 
Ino.14Gao.86As  layer  follows.  HRXRD  rocking  curves  were  obtained 
from  the  asymmetric  (444)  plane  using  the  to"  condition  at  +90°  and  - 
90°  specimen  rotation  angles  (see  chapter  2)  and  these  curves  are 
shown  in  Figure  5-1.  The  to"  peak  position  for  the  substrate  and 
epilayer  were  determined  from  the  two  curves  and  were  averaged  to 
obtain  the  true  to"  (co"  = 0 - <}>)  values  of  16.2156°  and  14.8361°  for 
substrate  and  epilayer,  respectively.  A similar  procedure  was 
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Figure  5-1.  X-ray  rocking  curves  obtained  using  the  co'  Bragg  condition 
at  (444)  reflections. 

a)  -90°  rotation  angle  and  b)  +90°  rotation  angle. 
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employed  in  the  determination  of  true  co+  (o)+  = 0 + $)  values  from  the 
scans  shown  in  Figure  5-2  which  were  found  to  be  125.5°  and 
124.3274°,  respectively,  for  substrate  and  epilayer.  The  equations  for 
a)-  and  co+  were  solved  to  obtain  the  values  of  0 and  <J>  for  substrate  and 
epilayer,  respectively.  These  values  were  used  in  the  equations  given 
in  Chapter  2 to  determine  the  perpendicular  and  in-plane  or  parallel 
lattice -mismatches  between  epilayer  and  substrate.  For  the  example 
shown  in  Figures  5-1  and  5-2,  the  perpendicular  and  parallel 
mismatches  were  found  to  be  10,200ppm  (1.02%)  and  6500ppm 
(0.65%),  respectively.  The  error  in  such  calculations  is  about  5%  of 
the  mismatch  value.  The  calculated  perpendicular  lattice-mismatch 
agreed  with  the  lattice-mismatch  determined  using  the  (004)  planes 
and  hence  the  assumption  of  tetragonal  distortion  appeared  to  be 
valid. 

The  theoretical  or  relaxed  lattice -mismatch  is  given  by 

epilayer  ~ ^GaAs  ^ (aInAs  ~ ‘klaAs^ 

V ^GaAs  'relaxed  ^GaAs 

where  x is  the  InxGai-xAs  composition  and  a is  the  lattice  parameter. 
The  residual  strain  in  the  layer  can  be  written  as 


17  Aa  ^ 

f Aa  1 1 

^GaAs 

.V^aAs  J|| 

^ 'relaxed  - 

3r 

where  the  parallel  symbol  denotes  the  in-plane  mismatch  that  is 
experimentally  determined  using  HRXRD  and  ar  is  the  fully  relaxed 
lattice  parameter  of  the  epilayer.  The  equation  is  normalized  to  the 
lattice  parameter  at  a particular  InxGai-xAs  composition.  For 
comparison  of  the  relaxation  between  single  layers  of  differing 
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Figure  5-2.  X-ray  rocking  curves  obtained  using  the  + Bragg  condition 
at  (444)  reflections, 
a)  -90°  rotation  and  b)  +90°  rotation. 
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compositions  a "percentage  relaxation"  of  the  layer  in  the  in-plane 
direction  is  defined  as 


% relaxation  = - 


1 - 


Aa 


■1- 


( Aa 


aGaAs  ;(|  V ^aAs  /relaxed  - 

f— 1 

V aGaAs  /relaxed 


x 100 


and  the  percent-relaxation  is  normalized  to  the  fully  relaxed  lattice- 
mismatch  strain.  In  order  to  obtain  the  relaxation  in  the 
perpendicular  direction,  the  perpendicular  lattice -mismatch  is  used 
instead  of  the  in-plane  mismatch  in  the  equation  given  above. 

Ino.14Gao.86As  layers  were  grown  to  various  thicknesses  and 
the  extents  of  relaxation  in  the  perpendicular  and  parallel  directions 
were  determined  as  a function  of  epilayer  thickness  (Figure  5-3).  The 
plot  in  Figure  5-3  shows  that  strain  relaxation  in  the  perpendicular 
direction  is  complete  at  a thickness  of  0.5pm.  The  strain  in  the 
perpendicular  direction  for  layer  thicknesses  smaller  than  0.5pm  was 
found  to  be  tensile  in  nature.  Relaxation  in  the  parallel  or  in-plane 
direction,  however,  was  not  complete  even  at  a thickness  of  2pm. 
Strain  relaxation  occurs  rapidly  up  to  a thickness  of  ~0.5pm,  as 
evidenced  by  the  steep  slope  in  the  percent  relaxation  versus 
thickness  curve  which  suggests  66%  strain  relief  at  0.5pm.  As  can  be 
seen  from  Figure  5-3,  relaxation  becomes  sluggish  with  further 
increase  in  layer  thickness,  showing  only  9%  further  relaxation  up  to  a 
thickness  of  2pm. 

InxGai-xAs  layers  having  compositions  in  the  range,  x=0.05-0.4, 
were  grown  to  a constant  thickness  of  1pm.  Also,  three  InxGai-xAs 
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Figure  5-3.  A plot  of  the  percent  relaxation  in  the  perpendicular  and 
parallel  directions  as  a function  of  thickness  in  a 
Ino.l4Gao.86As/GaAs  layer. 
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layers  with  x values  of  0.14,  0.27  and  0.4,  respectively,  were  grown  to 
a thickness  of  0.2pm  to  study  the  effect  of  composition  on  strain 
relaxation.  HRXRD  analysis  performed  on  these  layers  showed  that  the 
extent  of  relaxation  was  almost  a constant  (70%)  in  1pm  thick 
InxGai-xAs  layers  irrespective  of  composition  (Figure  5-4).  The 
thickness  at  which  the  relaxation  rate  (with  respect  to  thickness) 
becomes  slow  might  depend  on  the  misfit  strain  which  in  turn  is  a 
function  of  epilayer  composition.  Data  from  the  0.2pm  thick 
In^ai^As  layers  suggests  that  in  the  Ino.i4Gao.86As  layer  (small 
misfit  strain)  only  about  37%  of  the  mismatch  strain  is  relaxed 
whereas  in  Ino.27Gao.73As  and  Ino.4Gao.6As  layers  (large  misfit 
strains)  of  the  same  thickness  about  67%  of  the  strain  is  relaxed 
(Figure  5-4).  Thus  one  can  conclude  that  InxGai-xAs  layers  with 
larger  mismatch  strains  reach  the  plateau,  where  strain  relaxation  is 
s^u66ish.  at  much  smaller  thicknesses  compared  to  layers  with  smaller 
mismatch  strains.  Upon  reaching  this  plateau,  the  extent  of  strain 
relaxation  would  appear  to  be  a constant  for  any  layer  composition  as 
evidenced  by  the  data  shown  in  Figure  5-4. 

Currently  published  theories  on  strain  relaxation  [18,19,22] 
predict  values  of  residual  strain  as  a function  of  epilayer  thickness.  In 
order  to  compare  theory  with  experimental  data,  a plot  of  the  residual 
strain  versus  thickness  for  all  InxGai-xAs  layer  compositions  studied 
in  this  work  is  shown  in  Figure  5-5.  Strain  relaxation  theories  by 
Matthews  [3,4,18]  and  Tsao  et  al.  [19]  based  on  equilibrium  force 
balance  and  plastic  flow,  respectively,  suggest  that  residual  strain  is 
independent  of  InxGai-xAs  composition  and  only  a function  of  epilayer 


99 


Figure  5-4.  A plot  of  the  in-plane  relaxation  as  a function  of 

InxGai-xAs  composition,  x,  at  epilayer  thicknesses  of  0.2 
and  1pm. 


Residual  Strain  (ppm) 
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Thickness  (nm) 


Figure  5-5.  Residual  strain  data  obtained  for  various  InxGai-xAs 
compositions  and  plotted  as  a function  of  thickness 
showing  that  the  data  is  scattered  and  does  not  follow  any 
obvious  trend. 
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thickness.  However,  as  can  be  seen  from  Figure  5-5,  the  residual 
strain  is  a strong  function  of  InxGai-xAs  layer  composition,  ranging 
from  1000  to  8000ppm  with  increasing  composition  in  the  range 
x=0. 05-0.4.  Therefore,  it  can  be  concluded  that  these  composition 
independent  theories  are  invalid  in  the  InxGai-xAs/GaAs  material 
system. 

Maree  et  al.  [22]  proposed  a strain  relaxation  theory  based  on 
the  nucleation  of  surface  half-loops  and,  as  shown  in  Chapter  1 , 
experimental  data  from  various  studies  [9,22-24]  on  the 
In^ai^As/GaAs  system  appear  to  fit  the  theoretical  predictions 
within  an  error  of  about  2000ppm.  However,  for  the  composition 
range  studied  (x=0.07-0.25)  such  an  error  is  significant.  Even  though 
Maree  et  al.  [22]  suggest  that  their  theory  is  composition  independent, 
they  employ  a critical  half-loop  radius,  defined  as  that  radius  of  a half- 
loop which  maximizes  the  energy  of  the  loop,  which  is  a function  of 
the  critical  thickness  in  their  calculations.  The  critical  thickness,  in 
turn,  is  a function  of  InxGai-xAs  composition.  The  equation  for 
residual  strain  as  given  by  Maree  et  al.  [22]  is  as  follows: 


e = 


1 - v/2  1 ( 4rc  i \ 

1 + v h lln  b 2 J 


where  rc  is  the  critical  radius  of  the  half-loop  and  h is  the  thickness  of 
the  epilayer.  Maree  et  al.  [22]  state  that  the  thickness  of  the  layer 
must  be  maintained  such  that  h>rcsin<|)  where  <)>  is  the  angle  between 
the  glide-plane  of  a dislocation  and  the  plane  of  the  surface.  It  is  well 
known  that  strain  relaxation  occurs  when  the  thickness  of  the  epilayer 
exceeds  the  critical  thickness  (hc).  Therefore,  one  can  assume  that 
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for  the  surface  half-loop  nucleation  theory  to  hold  true,  hc=rcsin<f>. 
Since  the  critical  thickness  decreases  with  increasing  composition,  rc 
would  also  decrease  with  increasing  composition  since  <f>  remains  a 
constant  for  a specific  crystal  structure.  Such  a result,  however, 
suggests  that  residual  strain  in  the  InxGai-xAs  layer  would  decrease 
with  increasing  composition  (x)  which  is  contrary  to  the  experimental 
observations  reported  in  this  study  (Figure  5-5).  Thus,  Maree  et  al.'s 
[22]  surface  half-loop  nucleation  theory  cannot  explain  strain 
relaxation  in  the  InxGai-xAs/GaAs  material  system. 

5.3  Strain  Relaxation  in  Stacked  InvGai-vAs  Multilayer  Structures 

A schematic  of  a multilayer  structure  is  shown  in  Figure  5-6,  in 
order  to  facilitate  a description  of  the  results  obtained  in  this  case.  A 
stage  number  is  assigned  to  each  layer  forming  the  stacked  structure, 
the  numbers  increasing  with  increasing  distance  from  the  substrate. 
The  thickness  of  each  layer  was  0.2(im  unless  otherwise  specified. 
Separately  grown  multilayer  structures  comprising  n stages,  where 
n=l  to  8,  were  analysed  by  HRXRD  and  XTEM  with  the  results 
obtained  from  analysis  of  each  stage  being  used  as  the  basis  for  the 
design  of  each  subsequent  stage. 

Since,  as  discussed  in  the  above  section,  strain  relaxation  is  not 
complete  particularly  in  the  case  of  thin  InxGai-xAs  epilayers,  one  can 
think  of  an  unrelaxed  InxGai-xAs  layer  as  having  an  "effective 
composition"  which  will  always  be  less  than  the  actual  composition.  In 
this  work,  such  an  "effective  composition"  (xn')  for  the  InxnGai-xnAs 
layer,  is  defined  as  the  composition  corresponding  to  the  in-plane 
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Stage  8 (lnX8Gai_X8As) 
Stage  7 (lnX7Ga-|.X7As) 
Stage  6 (lnX6Ga-j.X6As) 
Stage  5 (In^Ga^xsAs) 
Stage  4 (lnX4Gai_X4As) 
Stage  3 (lnX3Ga-|.X3As) 
Stage  2 (lnX2Gai_X2As) 


Stage  1 (lnxiGa-j.x*jAs) 


Figure  5-6.  Schematic  of  the  stacked  InxGai-xAs  multilayer  system 

grown  by  molecular  beam  epitaxy  on  GaAs.  The  schematic 
shows  a total  of  eight  stages  (eight  InxGai-xAs  layers), 
however,  stacks  with  various  numbers  of  stages  were 
grown  and  analysed. 
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lattice  parameter  of  the  partially  relaxed  layer  (accounting  for  strain) 
which  is  evaluated  using  Vergard's  law  as  follows  : 

_ (an'  — ‘klaAs) 

(aInAs  — ^GaAs) 

Furthermore,  an  "effective  compositional  difference"  (Ax’)  between 
two  adjacent  stages  in  the  multilayer  structure,  for  example  stages  n 
and  (n+1),  is  defined  as  the  difference  between  the  relaxed 
composition,  x(n+l),  of  the  (n+l)**1  stage  and  the  effective 
composition,  xn',  of  the  n^1  stage. 

In  Chapter  3 the  notion  of  a "critical  composition"  (x=0.18)  for 
single  InxGai-xAs  layers  grown  on  GaAs  was  evoked.  In  this  section 
this  notion  is  extended  to  multilayer  structures  and  a "critical 
compositional  difference",  Ax=0. 18,  is  designated  with  respect  to  two 
adjacent  layers  in  a stack. 

Stage  1 of  the  multilayer  structure  comprised  an 
Ino.14Gao.86As,  i.e.  xl=0. 14,  layer  grown  on  GaAs.  Since  the  epilayer 
composition  is  below  the  "critical  composition"  (x=0.18),  the 
dislocations  are  confined  to  the  heterointerface  and  the  substrate  as 
previously  outlined  in  Chapter  3.  Using  HRXRD,  the  parallel  lattice- 
mismatch  between  Stage  1 and  the  GaAs  substrate  was  determined  to 
be  only  0.36%  which  is  well  below  the  theoretical  (completely 
relaxed)  lattice  parameter  mismatch  (Table  5-1)  of  1%,  as  evaluated 
assuming  Vergard's  law,  and  hence  Stage  1 is  only  partially  relaxed. 
The  experimentally  determined  lattice  parameter  corresponds  to  an 
"effective"  epilayer  composition  of  xl'=0.05,  compared  to  an  actual 
composition  of  0.14. 
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Table  5-1.  HRXRD  data  from  three  separately  grown  structures 

consisting  of  one-,  two-  and  three-stages,  respectively, 
denoting  the  experimentally  determined  in-plane  lattice- 
mismatch  of  the  appropriate  stages  in  each  structure  with 
respect  to  GaAs.  The  actual  compositions  of  the  three 
stages  were  xl=0.14,  x2=0.22  and  x3=0.3.  All  stages  were 
0.2pm  thick.  The  relaxed  lattice-mismatch  of  each  stage 
with  respect  to  GaAs  is  given  in  brackets  below  the 
corresponding  stage  number.  The  relaxed  compositions  of 
the  top  layers  are  denoted  in  brackets  next  to  their 
effective  compositions. 


Structure 

# 

Measured  in-plane 
lattice-mismatch  (%) 

Effective  composition 
of  the  top  layer 

Stage  1 Stage  2 Stage  3 

(1.0)  (1.58)  (2.15) 

1 

0.36 

0.05  (0.14) 

2 

0.80  1.0 

0.14  (0.22) 

3 

0.83  1.24  1.42 

0.20  (0.30) 
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Two  separate  structures  consisting  of  Inx2Gai-x2As  (Stage  2) 
compositions  of  x2=0.22  and  x2=0.27  were  grown  with  Stage  1 in 
each  case  being  x 1=0. 14  such  that  the  effective  compositional 
difference  between  the  top  layer  and  the  underlying  Inxl'Gai-xrAs 
(effective  composition)  layer  would  be,  in  one  case,  smaller  than  the 
critical  compositional  difference  (x2-xl'<0. 18)  and  in  the  other  case, 
greater  than  the  critical  compositional  difference  (x2-xl’>0.18). 

XTEM  analysis  of  the  structure  containing  the  Ino.22Gao.78As  (Stage 
2)  epilayer  showed  dislocation  evolution  at  the  Ino.22Gao.78As/ 
InO.14Gao.86As  (Stage2 /Stage  1)  heterointerface  and  in  the  underlying 
Ino.14Gao.86As  layer  (Figure  5-7a).  Also,  threading  dislocations  were 
not  observed  in  the  top  Ino.22Gao.78As  layer.  In  contrast,  XTEM 
analysis  of  the  structure  containing  the  Ino.27Gao.73As  (Stage  2) 
epilayer  revealed  threading  dislocations  in  the  top  Ino.27Gao.73As 
layer  as  well  as  dislocations  in  the  Ino.i4Gao.86As  (Stage  1) 
underlying  layer  (Figure  5- 7b). 

The  above  results  can  be  explained  as  follows.  If  the  effective 
compositional  difference  between  the  InxGai-xAs  layers  in  the  two 
stages  is  below  the  critical  compositional  difference  as  is  the  case  in 
the  structure  of  Figure  5.7a,  the  surface  image  force,  which  is  an 
inverse  function  of  layer  thickness,  would  not  be  significant  enough  to 
attract  dislocations  to  the  surface  due  to  the  larger  critical  thickness 
required  to  generate  misfit  dislocations.  Also,  if  the  yield  strength  of 
the  top  InX2Gai-x2As  layer  at  its  growth  temperature  is  greater  than 
that  of  the  underlying  InxiGai-xiAs  layer,  dislocation  evolution  would 
occur  at  the  two  heterointerfaces  and  in  the  underlying  InxiGai-xiAs 
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GaAs 


Figure  5-7.  Cross-sectional  transmission  electron  micrographs  of  two- 
stage  structures. 

a)  Effective  compositional  difference  smaller  than  the 
critical  compositional  difference  and  b)  Effective 
compositional  difference  greater  than  the  critical 
compositional  difference. 
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layer  and  the  GaAs  substrate  as  is  observed  in  the  micrograph. 

However,  if  the  effective  compositional  difference  between  the 
InxGai-xAs  layers  is  greater  than  the  critical  compositional  difference, 
as  is  the  case  in  the  structure  of  Figure  5-7b,  the  surface  image  force 
will  be  comparable  to  the  additional  forces  acting  on  a dislocation  at 
the  Inx2Gai-X2As/InxiGai-xiAs  heterointerface  since  the  critical 
thickness  in  this  case  will  be  small.  Thus,  the  dislocations  can 
traverse  to  the  surface  even  though  the  yield  strength  of  the  top  layer 
is  greater  than  that  of  the  underlying  layer.  The  composition  of 
x2=0.22  was  used  for  the  second  stage  in  all  subsequent  multilayer 
structures  since  no  threading  dislocations  were  observed  in  the 
structure  consisting  of  a Stage  2 having  this  composition. 

HRXRD  analysis  indicated  that  the  effective  composition  (taking 
into  account  partial  strain  relaxation)  of  the  selected  Stage  2 
In^Gai-x^As  layer  (x2=0.22)  is  x2'=0. 14  (Table  5-1).  Therefore,  an 
InX3Gai-x3As  composition  of  x3=0.3  was  chosen  for  Stage  3 of  the 
stacked  structure  since  x3-x2'<0.18.  By  the  same  reasoning,  a Stage  4 
composition  of  x4=0.35  was  selected  since,  as  shown  in  Table  5-1,  the 
effective  composition  of  Stage  3 (x3=0.3)  was  x3'=0.2  giving  again  an 
effective  compositional  difference  <0.18  at  the  Stage  3/Stage  4 
heterointerface.  XTEM  analysis  of  three-stage  and  four-stage 
structures  showed  a lack  of  threading  dislocations  in  the  top 
ImjcGai^As  layer  of  each  stack  as  illustrated  in  Figures  5-8a  and  5-8b. 
Strain  relief  occurs  by  the  formation  of  dislocations  at  the 
heterointerfaces  and  in  the  underlying  layers  (Stages  1-3).  Etch  pit 
analysis  of  the  Ino.35Gao.65As  top  layer  which  in  this  case  was  l(im 
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Figure  5.8:  Cross-sectional  transmission  electron  micrographs  of 

showing  dislocations  in  underlying  layers  and  a noticable 
absence  of  dislocations  in  the  top  InxGai-xAs  layer  of 
multilayer  structures. 

a)  three-stage  structure  and  b)  four-stage  structure. 
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thick  revealed  a threading  dislocation  density  of  only  l-2xl03/cm2  in 
this  material  (Figure  5-9). 

HRXRD  data  obtained  from  one-,  two-  and  three-stage  structures 
is  detailed  in  Table  5-1.  As  can  be  seen  from  the  table,  in  each  case 
the  top  layer  is  only  partially  relaxed.  Also,  strains  in  the  underlying 
InxGai-xAs  layers  approach  their  theoretical  values  with  increasing 
number  of  stages  suggesting  that  underlying  layers  continue  to 
undergo  strain  relief  by  the  generation  of  dislocations  upon  growth  of 
additional  stages.  It  is  also  noted  that  provided  the  thickness  of  the 
top  layer  and  the  effective  compositional  difference  with  respect  to 
the  underlying  layer  is  maintained  a constant,  the  extent  of  in-plane 
relaxation  of  the  top  InxGai-xAs  layer  with  respect  to  its  underlying 
layer  is  almost  a constant  as  evidenced  by  the  nearly  constant 
difference  between  the  actual  composition  and  the  effective 
composition  of  the  top  layer  following  each  stage  of  growth.  This 
observation  seems  to  suggest  a monotonic  increase  in  yield  strength  at 
the  growth  temperature  for  InxGai-xAs  with  increasing  values  of  x at 
least  up  to  x=0.3  since,  the  extent  of  relaxation  is  a direct  function  of 
yield  strength. 

Three  multilayer  structures  consisting  of  5,  6 and  7 stages, 
respectively,  were  grown  with  compositions  of  x5=0.4,  x6=0.45  and 
x7=0.48  above  Stage  4 (x4=0.35,0.2|im)  such  that  the  effective 
compositional  difference  between  the  top  two  stages  in  each  of  the 
three  structures  was  well  below  the  critical  compositional  difference. 
Stages  1-4  were  maintained  identical  in  each  case.  The  composition 
of  the  seventh  stage  was  chosen  to  be  x7=0.48  since  it  is  a 
composition  of  interest  for  potential  device  operation  at  1.55pm  and 
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Figure  5-9.  Scanning  electron  micrograph  showing  an  etch-pit  in  the 
top  Ino.35Gao.65As  layer  of  a four-stage  heterostructure. 
Dislocation  density  was  calculated  to  be  -lO^/cm^. 
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the  thickness  of  this  stage  was  selected  as  l|im  to  enable  etch  pit 
analysis  of  the  layer.  XTEM  analysis  of  the  multilayer  structures 
showed  that  strain  relief  continues  to  occur  by  the  generation  of 
dislocations  at  the  heterointerfaces  and  in  the  underlying  stages  as 
shown  in  Figures  5- 10a  and  5- 10b  which  show  the  structures 
containing  six  and  seven  stages,  respectively.  Etch  pit  analysis  of  the 
In0.48Gao.52As  layer  revealed  a threading  dislocation  density  of  only 
5-6xl03/cm2  (Figure  5-11).  These  observations  seem  to  indicate  that 
the  yield  strength  of  the  top  layer  at  the  growth  temperature  is  larger 
than  that  of  the  underlying  layer  at  each  stage  which  leads  to 
dislocation  formation  (yielding)  always  in  the  underlying  layers.  It  is  to 
be  noted  that  the  observed  dislocation  density  is  three  orders  of 
magnitude  smaller  than  that  obtained  by  using  the  multi-stage  strain 
relief  buffer  described  in  the  previous  chapter  even  though  the  top 
layer  composition  (x=0.48)  employed  in  this  heterostructure  is  larger 
than  that  in  the  top  layer  (x=0.4)  of  the  buffer  system  involving  SLSs. 

Finally,  a heterostructure  consisting  of  8 stages  was  grown  in 
which  stages  7 and  8 having  compositions,  x7=0.5  and  x8=0.53, 
respectively,  were  grown  above  Stage  6 (x6=0.45).  Stages  1-5  were 
identical  to  those  involved  in  the  structure  consisting  of  five  stages. 
XTEM  analysis  of  the  eight  stage  structure  revealed  threading 
dislocations  in  the  Ino.53Gao.47As  top  layer  as  shown  in  Figures  5- 12a 
and  5- 12b.  Dislocations,  however,  were  not  observed  in  the 
underlying  Ino.5Gao.5As  (Stage  7)  layer.  Since  the  effective 
compositional  difference  between  the  top  two  layers  is  below  the 
critical  compositional  difference,  the  presence  of  dislocations  in  the 
top  layer  and  not  its  underlying  layer  is  explained  if  the  yield  strength 
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Figure  5-10.  Cross-sectional  transmission  electron  micrographs  of 
multilayer  structures  showing  a lack  of  threading 
dislocations  in  the  top  layer  of  each  structure  and 
extensive  dislocation  evolution  in  the  underlying  layers, 
a)  six-stage  structure  and  b)  seven-stage  structure. 
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Figure  5-1 


♦Vi {^4 


Scanning  electron  micrograph  showing  etch-pits  in  the 
top  Ino.48Gao.52As  layer  of  a seven-stage  structure. 
Dislocation  density  was  calculated  to  be  5-6xl()3/cm2  in 
this  material 
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Figure  5-12.  Cross-sectional  transmission  electron  micrographs  of  an 
eight-stage  structure. 

a)  Showing  almost  the  whole  structure  with  the 
exception  of  the  GaAs  substrate  and  the  Ino.i4Gao.86As 
layer; 
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Figure  5-12— continued 

b)  An  Ino.5Gao.5As  (Stage  7)  layer  sandwiched  between 
the  Ino.53Gao.47As  (Stage  8)  and  Ino.45Gao.55As  (Stage 
6)  layers.  The  micrograph  shows  an  apparant  lack  of 
dislocations  in  the  Ino.5Gao.5As  layer  and  threading 
dislocations  in  the  Ino.53Gao.47As  layer.  The  critical 
compositional  difference  was  not  exceeded  at  any  of  the 
heterointerfaces,  however,  the  yield  strength  of 
Ino.53Gao.47As  at  the  growth  temperature  would  appear 
to  be  smaller  than  that  of  Ino.5Gao.5As,  resulting  in 
dislocation  evolution  in  the  Ino.53Gao.47As  material. 
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at  the  growth  temperature  of  Ino.53Gao.47As  is  smaller  than  that  of 
Ino.5Gao.5As.  The  absence  of  dislocations  in  the  Ino.5Gao.5As  layer 
suggests  that  its  yield  strength  is  greater  than  that  of  the  underlying 
Ino.45Gao.55As  layer. 

Studies  reported  in  the  literature  [65,661  on  the  high 
temperature  mechanical  properties  of  InAs  and  GaAs  using 
compression  tests  have  shown  that  InAs  has  a lower  yield  strength 
than  GaAs  at  an  elevated  temperature  (540°C).  Combining  the 
dislocation  evolution  information  obtained  in  this  work  and  available 
data  on  the  high  temperature  yield  strengths  for  binary  InAs  and  GaAs, 
one  can  sketch  the  yield  strength  (in  arbitrary  units)  of  InxGai-xAs  as 
a function  of  In  content  with  respect  to  the  particular  temperatures 
employed  for  alloy  growth  as  shown  in  Figure  5-13.  For  the  growth 
temperatures  employed  in  this  study,  a peak  in  the  yield  strength 
would  seem  to  occur  at  an  In  content  of  x=0.5.  Therefore,  for  x<0.5, 
one  can  use  the  notion  of  "critical  compositional  difference"  to  obtain 
low-dislocation  density  In^ai-jcAs  alloys  grown  on  GaAs,  since  strain 
relief  occurs  by  the  generation  of  dislocations  in  the  softer  (lower  yield 
strength)  material  i.e.  the  underlying  layer  in  the  structure,  and,  since 
below  the  critical  compositional  difference  the  surface  image  force  is 
not  large  enough  to  influence  the  dislocations  at  the  heterointerface. 

In  this  study,  all  of  the  InxGai-xAs  sacrificial  layers  in  the 
multilayer  system  were  grown  to  a thickness  of  0.2pm.  However, 
based  on  the  residual  strain  studies  performed  on  single  InxGai-xAs 
layers  grown  on  GaAs  detailed  in  the  previous  section,  it  should  be 
possible  to  employ  fewer  stages  in  the  multi-stage  system  in  order  to 
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Figure  5-13.  A conceptual  plot  of  the  yield  strength  at  the  growth 
temperature  versus  InxGai-xAs  composition  based  on 
indirect  evidence  from  dislocation  evolution  and  strain 
relief  studies  (see  text).  The  InxGai-xAs  layers  employed 
in  this  study  and  their  corresponding  growth 
temperatures  are  as  follows:  x=0.14,  520°C;  x=0.22, 
450°C;  x=0.3,  420°C;  x=0.35,  400°C;  x=0.4,  380°C; 
x=0.45,  360°C;  x=0.48,  350°C:  x=0.5.  340°C  and  x=0.53. 
330°C. 
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achieve  a particular  InxGai-xAs  composition  by  growing  thicker 
sacrificial  layers  than  were  grown  in  this  study. 


CHAPTER  6 

CONCLUSIONS  AND  RECOMMENDATIONS 

Several  conclusions  can  be  drawn  from  the  work  described  in 
this  dissertation  which  are  as  follows: 

a)  A "critical  composition"  of  x=0. 18  exists  for  the  InxGai-xAs/  GaAs 
heteroepitaxial  system  above  which  threading  dislocations  propagate 
into  the  InxGai-xAs  epilayer  and,  in  general,  dislocation  evolution  in 
the  InxGai-xAs/GaAs  system  as  a function  of  composition  can  be 
explained  employing  a force  balance  model  which  includes  a surface 
image  force  term, 

b)  InxGai-xAs  layers  with  relatively  low-dislocation  densities 
(=10® /cm2)  can  be  grown  indirectly  on  GaAs  up  to  a composition  of 
x=0.4  using  multistage  buffer  structures  consisting  of  conventional 
dislocation  blockers,  i.e.  strained -layer  superlattices.  The  inability  of 
such  multistage  buffer  structures  to  reduce  the  dislocation  density  in 
InxGai-xAs  layers  having  x>0.4  appears  to  be  due  to  dislocation 
nucleation  and  evolution  within  the  SLSs  comprising  the  buffer. 

c)  Strain  relaxation  in  InxGai-xAs/GaAs  layers  is  complete  in  the 
perpendicular  (growth)  direction  at  relatively  small  layer  thicknesses 
(~0.5pm),  however,  a finite  (-25%)  residual  strain  exists  in  the  in- 
plane direction  even  for  relatively  thick  layers.  Residual  strain  in  the 
InxGai-xAs  layers  is  strongly  dependent  on  both  the  layer  thickness 


120 


121 


and  composition,  decreasing  with  increasing  thickness  and  increasing 
with  increasing  composition. 

d)  Extremely  low-dislocation  density  (<104/cm2)  InxGai-xAs  layers 
can  be  grown  up  to  an  x value  of  ~0.5  indirectly  on  GaAs  by  use  of 
multilayer  epitaxial  structures  that  do  not  involve  conventional  buffers. 
Such  structures  are  designed  such  that  the  "effective  compositional 
difference"  (accounting  for  partial  relaxation)  between  successive 
In^ai-jcAs  layers  in  the  stack  is  below  a "critical  compositional 
difference"  of  x=0.18. 

e)  A monotonic  increase  in  the  elevated  temperature  yield  strength  of 
InxGai-xAs  as  a function  of  composition  is  suggested,  based  on 
indirect  evidence,  for  x values  in  the  range,  0 to  0.5,  while  the  alloy 
yield  strength  is  suggested  to  decrease  with  increasing  composition 
for  x values  >0.5  for  the  growth  temperatures  employed  in  this  study. 

The  above  conclusions  suggest  several  possibilities  for  future 
research,  both  fundamental  and  applied,  in  the  field  of  compound 
semiconductor  heteroepitaxy.  On  a fundamental  note,  the  present 
work  has  highlighted  the  need  for  data  on  the  elevated  temperature 
mechanical  properties  of  compound  semiconductors,  particularly 
ternary  and  quartemary  alloys.  Such  data  could  indicate  the  possibility 
of  growing  device  quality  alloys  in  other  lattice -mismatched  systems 
assuming  that  they  would  behave  in  a manner  similar  to  that  displayed 
by  the  InxGai-xAs/GaAs  system,  i.e.,  the  alloy  would  have  to  be  harder 
at  the  appropriate  growth  temperature. 

On  the  applied  research  front,  dislocation-sensitive  devices  such 
as  light  emitting  diodes  and  diode  lasers  could  be  designed  based  on 
large  In  content,  low-dislocation  density  InxGai-xAs  material.  The 
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testing  of  such  devices  should  further  demonstrate  the  high-quality 
nature  of  the  InxGai-xAs  material  which  has  been  achieved  as  a result 
of  the  fundamental  work  described  in  this  dissertation.  A more  far 
reaching  research  area  would  be  the  design  of  monolithically 
integrated  photonic  devices  operating  at  1.3  and  1.55pm  with  GaAs 
ICs. 
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